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Abstract;  Within  the  System-Wide  Water  Resources  Program  (SWWRP), 
multiple  riverine,  estuarine,  watershed,  and  subsurface  flow  models  are 
being  modified  to  address  issues  of  environmental  concern.  Several 
integration  approaches  are  either  ongoing  or  proposed  to  accomplish  this 
task.  To  have  a  full  system-wide  water  quality  and  contaminant  capability 
in  SWWRP,  the  different  hydrologic  and  hydraulic  engines  must  utilize  a 
common  water  quality  and  contaminant  approach  to  prevent  the  arbitrary 
partitioning  of  constituents.  The  goal  of  this  development  effort  has  been 
to  upgrade  existing  hydrologic  and  hydraulic  models  (i.e.,  water  engines) 
using  a  common  water  quality  approach  in  order  to  facilitate  their  linkage 
and  application  on  a  system-wide  basis. 

In  keeping  with  a  common  water  quality  approach  to  model  development, 
a  library  of  water  quality  kinetics  has  been  developed  such  that  these 
kinetics  can  be  integrated  with  a  variety  of  water  transport  engines.  The 
library  of  algorithms  are  able  to  deal  with  a  multi-species,  multi-phase, 
multi-reaction  system,  should  include  both  fast  (equilibrium-based)  and 
slow  (non-equilibrium-based  or  rate-based)  reactions,  are  easily  extensible 
to  new  reaction  pathways,  include  both  common  nutrient  and 
contaminant  packages  as  well  as  geochemistry,  and  have  a  simple,  well- 
defined  data  interface  and  calling  procedure,  making  them  portable.  The 
modules  are  developed  such  that  they  are  data  structure  independent  thus 
facilitating  their  integration  into  a  wide  range  of  modeling  systems. 
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Preface 

This  report  summarizes  the  development  of  overland  and  channel  nutrient 
kinetic  modules  for  linkage  with  a  variety  of  hydraulic  and  hydrologic 
modeling  systems.  This  development  effort  was  performed  by  the  Engi¬ 
neer  Research  and  Development  Center  (ERDC),  Waterways  Experiment 
Station,  Vicksburg,  Mississippi.  Eunding  was  provided  under  the  System- 
Wide  Water  Resources  Program  (SWWRP).  Dr.  Steven  L.  Ashby  is  the 
program  manager  for  SWWRP.  Appreciation  is  extended  to  all  who 
assisted  in  the  formulations  and  review  of  the  kinetics  implemented  within 
the  nutrient  sub-modules  (NSM). 

Principal  Investigators  for  this  study  were  Dr.  Billy  E.  Johnson  of  the 
Water  Quality  and  Contaminant  Modeling  Branch  (WQCMB)  and  Terry  K. 
Gerald  of  SpecPro  Inc.  Mr.  Gerald  was  funded  as  an  onsite  contractor 
under  Task  Order  Contract  W912HZ-05-D-0011  on  Civil  Delivery  Order 
No.  15.  Dr.  Johnson  conducted  his  portion  of  the  study  under  the  general 
supervision  of  Dr.  Barry  W.  Bunch,  Chief  of  WQCMB,  and  under  the  gen¬ 
eral  supervision  of  Dr.  Richard  E.  Price,  Chief  EPED,  and  Dr.  Beth  C. 
Eleming,  Director  of  EL. 

COL.  Richard  B.  Jenkins  was  Commander  and  Executive  Director  of 
ERDC.  Dr.  James  R.  Houston  was  Director. 
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Summary 

Nutrient  Sub-Model  (NSM)  version  i.o  consists  of  Overland  and  Channel 
Nutrient  Kinetics.  These  kinetics  were  initially  taken  from  SWAT 
formulations;  however,  modifications  were  made  in  an  effort  to  enhance 
their  linkage  with  USAGE  hydraulic  and  hydrologic  water  engines. 

The  overland  module  simulates  the  nitrogen  and  phosphorus  nutrient 
cycles  with  algorithms  available  for  initializing  soil  concentrations  if 
observed  data  are  not  available.  The  channel  module  simulates  nitrogen 
and  phosphorus  in  addition  to  accounting  for  algal  growth,  carbonaceous 
biological  oxygen  demand,  and  dissolved  oxygen. 

Initial  integration  is  currently  taking  place  with  the  Gridded  Surface  Sub¬ 
surface  Hydrologic  Analysis  (GSSHA)  model  and  the  Hydrologic  Engineer 
Center  River  Analysis  System  (HEC-RAS)  model.  GSSHA  is  a  physically 
based  distributed  watershed  model  that  will  take  advantage  of  all  the  NSM 
capabilities.  HEC-RAS  is  a  i-D  stream  model  that  will  only  use  the  channel 
kinetic  capabilities  of  NSM. 

Euture  versions  of  NSM  will  include:  i)  a  carbon  cycle;  2)  plant/soil 
dynamics;  and  3)  updated  channel  kinetics.  As  these  new  features  are 
implemented,  an  updated  reference  manual  will  be  created  and  integration 
into  various  H&H  modeling  systems  will  be  done.  In  addition,  as  field 
studies  and  new  understandings  of  the  nutrient  and  water  quality  proc¬ 
esses  are  performed  and  discovered,  modifications  to  the  NSM  algorithms 
will  take  place  and  be  propagated  to  the  appropriate  modeling  systems. 

The  common  framework  approach  taken  in  this  development  effort  will 
ensure  that  the  latest  nutrient  process  descriptions  are  implemented  into 
the  USAGE  modeling  systems  in  the  most  efficient  manner. 
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1  Introduction 

Background 

The  concept  of  watershed  planning  is  not  new  to  the  U.S.  Army  Corps  of 
Engineers.  Throughout  its  history,  the  Corps  has  incorporated  watershed 
planning  into  the  process  by  which  it  manages  water  resource  systems. 
Even  the  Corps’  geographic  organization,  along  watershed  boundaries 
rather  than  State  and  county  lines  in  most  cases,  supports  the  historic 
understanding  of  the  need  to  manage  water  within  a  watershed  context. 
However,  this  understanding  and  organizational  concept  alone  are  not  suf¬ 
ficient  to  ensure  proper  protection  and  responsible  development  of  the 
Nation's  water  resources  in  the  future. 

This  country  is  facing  a  looming  water  crisis.  It  is  seeing  frequent  regional 
droughts,  disputes  over  allocation  brought  on  by  growing  population 
demands,  environmental  degradation  due  to  changing  land  uses,  and 
widespread  disagreement  over  competing  purposes  for  water  resources 
use.  The  problems  of  rapid  growth  in  certain  areas  are  worse  because 
responsibilities  to  address  water  needs  are  distributed  among  a  multitude 
of  government  agencies  and  private  companies,  so  the  problem-solving 
efforts  are  typically  fragmented.  The  results  are  predictable:  instead  of 
broadly  supported  regional  solutions  that  address  multiple  needs,  balance 
competing  uses,  and  can  be  quickly  implemented,  solutions  are  narrowly 
focused,  contentious  and  slowly  implemented,  uncertain,  expensive,  and 
inferior.  Eor  several  years,  the  Nation's  priorities  and  values  related  to 
water  resources  have  been  changing.  This  is  a  natural  evolution  resulting 
from  advances  in  scientific  knowledge,  public  reaction  to  that  increase  in 
knowledge,  and  an  unprecedented  national  prosperity  that  allows  us  to 
consider  more  than  the  short-term  basics  of  life. 

In  the  summer  and  fall  of  2000,  the  Corps  of  Engineers  held  a  series  of  16 
“listening  sessions”  around  the  Nation  to  hear  what  Americans  thought 
were  the  major  water  challenges  for  the  21st  Century.  The  participants 
provided  valuable  input  for  Eederal  involvement  that  would  best  help 
various  levels  of  government  face  these  challenges.  One  of  the  frequently 
raised  topics  was  the  need  to  address  water  challenges  from  a  watershed 
view,  highlighting  collaboration  and  integration.  Some  present-day 
watershed  management  efforts,  such  as  the  Comprehensive  Everglades 
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Restoration  Plan,  already  promote  active  participation  of  all  interested 
parties  in  the  planning  and  decision-making  process.  The  Corps  believes 
that  this  concept  of  integration  is  the  key  to  reforming  America's  water 
development,  protection,  and  restoration.  In  its  recently  released 
Watershed  Perspeetivefor  the  Civil  Works  Program,  the  Corps  describes 
the  foundation  for  watershed  activities  and  involvement.  The  nine 
Watershed  Principles  outlined  there  provide  the  approach  the  Corps  seeks 
to  follow  in  its  water  resource  management. 

•  Seeking  sustainable  water  resources  management. 

•  Integrating  water  and  related  land  management. 

•  Considering  future  water  demands. 

•  Coordinating  planning  and  management. 

•  Promoting  cooperation  among  government  agencies  at  all  levels. 

•  Encouraging  public  participation. 

•  Evaluating  monetary  and  non-monetary  tradeoffs. 

•  Establishing  interdisciplinary  teams. 

•  Applying  adaptive  management  as  changing  conditions  or  objectives 
warrant. 

Unlike  the  single-purpose,  project-driven  initiatives  that  the  Corps  has 
been  directed  to  accomplish  in  the  past,  the  perspective  of  this  new  water¬ 
shed  approach  is  based  on  multi-purpose,  multi-objective  management, 
examining  all  water  needs  in  the  watershed  and  receiving  water  bodies. 
With  this  broader  context,  watershed  partners  would  collaborate  to  simul¬ 
taneously  address  multiple  objectives  -  environmental  quality,  social 
effects,  and  national  and  regional  economic  development. 

In  support  of  the  Corps’  watershed  approach,  the  System- Wide  Water 
Resources  Program  (SWWRP)  was  designed  to  assemble  and  integrate  the 
diverse  components  of  water  resources  management.  Products  from  this 
program  are  designed  to  help  users  surpass  individual  project  level 
analysis,  and  apply  current  and  improved  technologies  for  multi¬ 
disciplinary  system-wide  assessments.  The  ultimate  goal  of  SWWRP  is  to 
provide  the  Corps,  its  partners,  and  stakeholders  with  the  overall  techno¬ 
logical  framework  and  analytical  tools  to  restore  and  manage  water 
resources  and  balance  human  development  activities  with  natural  system 
requirements. 
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Within  SWWRP,  multiple  riverine,  estuarine,  watershed,  and  subsurface 
flow  models  are  being  modified  to  address  issues  of  environmental  con¬ 
cern.  Several  integration  approaches  are  either  ongoing  or  proposed  to 
accomplish  this  task.  To  have  a  full  system-wide  water  quality  and  con¬ 
taminant  capability  in  SWWRP,  the  various  hydrologic  and  hydraulic 
engines  must  utilize  a  common  water  quality  and  contaminant  approach  to 
prevent  the  arbitrary  partitioning  of  constituents.  The  goal  of  this  develop¬ 
ment  effort  is  to  upgrade  existing  hydrologic  and  hydraulic  models  (i.e., 
water  engines)  using  a  common  water  quality  approach  in  order  to  facili¬ 
tate  their  linkage  and  application  on  a  system-wide  basis. 

In  keeping  with  a  common  water  quality  approach  to  model  development, 
a  library  of  water  quality  kinetics  will  be  developed  such  that  they  can  be 
integrated  with  a  variety  of  water  transport  engines.  The  library  of  algo¬ 
rithms  will  be  able  to  deal  with  a  multi-species,  multi-phase,  multi¬ 
reaction  system,  should  include  both  fast  (equilibrium-based)  and  slow 
(non-equilibrium-based  or  rate-based)  reactions,  will  be  easily  extensible 
to  new  reaction  pathways,  will  include  both  common  nutrient  and  con¬ 
taminant  packages  as  well  as  geochemistry,  and  will  have  a  simple,  well- 
defined  data  interface  and  calling  procedure,  making  them  portable.  The 
modules  will  be  developed  such  that  they  are  data  structure  independent, 
thus  facilitating  their  integration  into  a  wide  range  of  modeling  systems. 

Objectives 

The  primary  purpose  of  the  System- Wide  Water  Resources  Program 
Nutrient  Sub-Model  (SWWRP-NSM)  is  to  provide  a  method  of  deter¬ 
mining  the  fate  of  nitrogen,  phosphorus,  and  carbon  in  watersheds 
exposed  to  receiving  water.  This  will  be  integrated  within  the  distributed 
watershed  model’s  solute  transport  component,  so  the  concentrations  of 
the  nutrient  species  are  simulated  for  every  watershed  grid,  and  are 
updated,  along  with  water  flow  and  solute  transport,  every  time-step. 

This  report  describes  the  current  state  of  the  NSM  in  addition  to  the  nutri¬ 
ent  process  theory  and  formulations  proposed  for  future  versions  of  NSM. 
The  NSM  simulates  three  domain  processes:  soil,  overland  flow,  and 
stream/channel.  Chapter  2  discusses  the  soil  nutrient  cycle  processes. 
Chapter  3  discusses  the  overland  flow  nutrient  transport  and  transforma¬ 
tion  processes,  and  Chapter  4  discusses  the  in-stream  nutrient  transport 
and  water  quality  processes. 
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As  research  continues,  it  is  anticipated  that  improved  process  descriptions 
will  be  developed  and  as  such  will  be  integrated  into  the  NSM.  As  these 
development  efforts  occur,  this  report  will  be  revised  in  order  to  reflect  the 
current  options  available  within  the  NSM. 
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2  Soil  Nutrient  Processes 

Nitrogen  (N)  and  phosphorus  (P)  are  essential  nutrients  for  plants  and 
animals;  however,  excessive  accumulation  of  nutrients  can  represent  too 
much  of  a  good  thing  and  cause  eutrophication,  especially  in  water  bodies 
such  as  streams,  lakes,  and  estuaries.  The  problem  of  eutrophication  has 
drawn  attention  to  the  role  of  watershed  management  in  elevating  N  and  P 
concentrations  within  soil  leachates  and  runoff.  The  main  factors  influenc¬ 
ing  watershed  nutrient  movement  can  be  separated  into  the  transport, 
nutrient  source,  and  management  factors.  Transport  factors  include  the 
mechanism  by  which  nutrients  move  in  runoff  and  erosion  within  the 
landscape.  Factors  which  influence  the  source  and  amount  of  nutrients 
available  to  be  transported  are  soil  nutrient  content  and  form  of  nutrient 
applied.  Nutrient  management  factors  include  the  method  of  application, 
timing,  and  placement  in  the  landscape  as  influenced  by  the  management 
of  application  equipment  and  tillage.  Effective  nutrient  management, 
whether  directed  toward  nutrient  supply  or  abatement,  requires  a  working 
knowledge  of  N  and  P  cycling  that  involves  a  series  of  transformation  and 
transport  processes. 

The  complexity  of  nutrient  transformations  and  transport  in  watershed 
systems  has  led  to  the  use  of  both  conceptual  and  mechanistic  models  for 
their  description.  Modeling  of  nutrients  in  the  watershed  consists  of  three 
distinct  parts.  The  first  part  deals  with  simulating  most  of  the  nutrient 
transformations  and  movements  in  the  soil  nutrient  cycle,  whereas  the 
second  part  focuses  on  the  transport  and  transformation  of  nutrients  in 
the  overland  flow.  The  third  part  simulates  nutrient  in-stream  transport 
and  transformation  processes.  The  distribution  of  soil  nutrients  within  a 
watershed  can  provide  clues  as  to  how  nutrients  are  transported  from  soil 
to  stream.  The  Nutrient  Sub-Model  (NSM)  has  been  developed  such  that  it 
simulates  the  nutrient  processes,  in  soil,  overland  flow  and  stream,  for 
multiple  species  and  phases. 

Within  NSM,  soil  nutrient  cycles  will  be  modeled  for  N  and  P.  Mathemati¬ 
cal  models  of  the  soil  N  and  P  are  generally  in  the  form  of  storage  (pool) 
accounting  procedures;  however,  in  a  very  limited  number  of  models  they 
have  been  modeled  by  a  set  of  partial  differential  equations  derived  from 
the  actual  physical  and  chemical  laws  which  describe  the  transformation 
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process.  This  pool  conceptual  model  of  N  and  P  in  soils  is  relatively 
straightforward.  In  comparison  to  other  nutrient  processes  models,  the 
soil  N  and  P  representation  in  the  SWAT  model  (Neitsch  et  al.  2002)  offers 
a  fair  accounting  for  most  of  the  variables  in  the  soil  N  and  P  cycle  and  it 
can  provide  an  adequate  estimate  of  N  and  P  loads  transported  with  the 
overland  flow  in  dissolved  or  particulate  forms. 

Although  the  dynamics  of  the  soil  N  and  P  cycle  have  been  described  ade¬ 
quately  in  the  SWAT  model,  the  relations  in  the  model  have  not  been  for¬ 
mulated  in  partial  differential  form.  Nevertheless  these  relations  can  be 
utilized  to  build  a  more  mathematically  sound  model  to  simulate  the  soil  N 
and  P  dynamics  in  a  partial  differential  equation  formulation.  In  NSM  the 
soil  nutrient  mass  balance  equations  for  each  state  variable  define  all  the 
inputs,  outputs  or  storage  of  nutrients  in  a  system,  or  in  a  specified  region. 
The  mathematical  formulations  describing  the  transformation  and  the  bio¬ 
chemical  processes  of  the  different  soil  N  and  P  variables  in  SWAT  were 
adopted  along  with  recent  research  findings  to  develop  NSM.  The  work 
represented  here  aims  to  improve  the  spatial  resolution  and  the  level  of 
physical  description  in  nutrient  modeling  by  the  development  of  a  grid- 
based  nutrient  modeling  component  to  allow  the  modeling  of  nutrients 
using  fully  distributed  watershed  models.  All  nutrient  mass  balance  equa¬ 
tions  are  based  on  watershed  cells  and  maintained  for  multiple  soil  layers. 
Most  of  the  possible  chemical  and  physical  processes  occurring  to  nutri¬ 
ents  in  both  the  soil  and  the  overland  flow  have  been  accounted  for  in  the 
NSM.  The  hydrological  variables  required  to  drive  the  NSM  can  be  calcu¬ 
lated  using  any  physically  based  distributed  model  capable  of  producing  a 
reasonable  simulation  of  the  flow  and  sediment  fields  in  the  watershed. 
Linking  a  physically  based  hydrological  model  with  NSM  gives  a  more 
realistic  description  of  the  nutrient  dynamics  in  watersheds.  This  is  espe¬ 
cially  important  for  agricultural  watersheds  where  the  nutrient  plays  a 
more  important  role  and  its  occurrence  is  highly  variable  both  in  time  and 
space.  Hence,  nutrient  contamination  hot-spots  are  more  accurately  iden¬ 
tified  and  watershed  management  changes  to  reduce  nutrient  transport 
can  be  made  more  confidently. 

In  order  to  better  describe  the  interaction  between  mineral  and  organic 
N/P  pools,  NSM  will  include  a  module  that  simulates  carbon  (C)  dynamics 
in  the  soil-plant  system  according  to  a  C  cycle  similar  to  that  used  in  the 
EDYS  model  (Childress  et  al.  2002).  It  includes  three  soil  organic  matter 
(SOM)  pools:  active,  slow,  and  passive.  The  carbon  cycle  is  fundamental 
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for  simulation  of  all  organic  matter  dynamics  and  many  nutrient  cycling 
processes.  To  be  effective  in  real-world  applications,  NSM  also  must  be 
connected  with  models  for  plant  growth,  soil  temperature,  water  and  sol¬ 
ute  transport,  as  well  as  other  factors. 

Soil  nitrogen  module 

Nitrogen  cycle 

Of  all  the  mineral  nutrients,  N  has  the  most  complex  nutrient  cycle,  largely 
because  N  can  exist  as  a  gas  (both  ammonia  and  nitrogen  gas),  whereas 
the  other  mineral  nutrients  do  not  exist  as  gases  under  normal  soil  condi¬ 
tions.  To  help  understand  the  various  components  of  the  N  cycle,  defini¬ 
tions  and  molecular  formulas  of  the  numerous  N  forms  are  provided  in 
Table  i. 


Table  1.  Definitions  of  each  N  form. 


Nitrogen 

Form 

Molecular 

Formula 

Notes 

Nitrogen  gas 

N2  (g) 

- 

Ammonia  gas 

NH3  (g) 

- 

Ammonium 

NHl 

Plant  available,  attracted  to  exchange  sites  on  clay 
particles 

Nitrate 

NO3 

Very  mobile,  requires  more  energy  by  plant  than 
ammonium 

Nitrite 

NO2 

Mobile,  generally  low  concentrations 

Organic  N 

- 

Slowly  supplies  available  N  to  soil  solution 

N  occurs  in  several  chemical  forms.  These  forms  vary  greatly  not  only  in 
their  characteristics  but  also  in  the  way  they  behave.  These  forms  are  clas¬ 
sified  as  either  organic  or  inorganic  N.  Typically,  most  N  in  soils  and  surfi- 
cial  sediments  occurs  in  organic  form.  Organic  forms  of  N  are  found  in 
compounds  such  as  amino  acids,  protein,  and  more  resistant  N  com¬ 
pounds  (ultimately,  humus).  Inorganic  forms  include  ammonium  (NH^), 

nitrate  ( NO~ ),  and  with  a  low  concentration  of  nitrite  ( NO^ ).  Organic  N 
and  NH'^  are  mostly  absorbed  by  clays.  In  such  forms,  it  can  be  consid¬ 
ered  immobile.  Inorganic  forms  of  N  are  the  “available”  forms  that  plants 
and  microorganisms  can  use  or  that  can  move  in  the  soil  as  water  moves 
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through  it.  Most  N  in  the  soil  is  in  the  unavailable  organic  form.  Many 
forms  of  N  are  present  in  any  environment  at  any  point  in  time,  mainly 
because  N  readily  shifts  from  one  form  to  another.  NO~  and  NH'^  can  be 

converted  to  organic  N  (ON)  by  plants  and  microorganisms;  ON  can  be 
converted  back  to  inorganic  forms  as  the  organic  compounds  decompose. 
N  can  also  shift  between  inorganic  forms.  These  shifts  occur  as  nature 
attempts  to  establish  an  equilibrium  among  the  various  forms  as  environ¬ 
mental  conditions  change.  N  in  the  soil  interacts  with  the  atmosphere,  soil 
particles,  soil  solution,  microorganisms,  and  plants.  If  a  new  source  of  N  is 
added  to  alter  the  balance  or  if  environmental  conditions  (such  as 
temperature  and  moisture)  change,  N  transformations  take  place.  Because 
environmental  conditions  are  constantly  changing,  N  transformations  are 
constantly  occurring.  This  continual  movement  of  N  from  one  form  to 
another  is  known  as  the  “nitrogen  cycle.”  The  N  cycle  consists  of  various 
storage  pools  of  N  and  processes  by  which  the  pools  exchange  N.  Figure  i 
shows  the  major  components  of  the  watershed  model  N  cycle. 


Soil  N  cycling  consists  of  nine  major  processes:  mineralization,  immobili¬ 
zation,  nitrification,  denitrification,  volatilization,  sorption,  plant  uptake. 
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leaching,  and  surface  runoff  and  erosion.  Each  of  these  processes  is 
described  below. 

Soil  nitrogen  transformations 

Nitrogen  exists  in  a  number  of  chemical  forms  and  undergoes  chemical 
and  biological  reactions.  Unlike  other  nutrients,  only  a  small  amount  of  N 
is  contributed  from  the  mineral  part  of  the  soil  (i.e.,  rocks).  Most  naturally 
occurring  N  enters  the  soil  either  as  NO~  or  NH'^  in  rainfall  (atmospheric 

deposition)  or  by  special  plants,  such  as  alfalfa,  ceanothus,  and  red  alder, 
that  are  called  “nitrogen  fixers.”  Humans  also  increase  the  N  in  soils  by 
fertilizing  with  either  chemical  or  organic  fertilizers.  Once  in  the  soil,  the  N 
will  transform  through  the  processes  of  mineralization,  immobilization, 
volatilization,  nitrification,  denitrification,  plant  uptake,  and  sorption. 

Two  of  these  transformations— volatilization  and  denitrification— result  in 
losses  of  N  from  the  soil.  N  is  also  lost  through  leaching  and  runoff. 
Because  most  of  the  transformations  and  losses  are  caused  by  microbial 
action,  all  the  processes  are  slowed  down  considerably  when  temperatures 
are  low. 

Mineralization  and  immobilization 

Mineralization  of  N  occurs  when  the  organic  matter  (OM)  in  biosolids 
decomposes.  The  soil  microorganisms  break  organic  bonds  to  obtain 
energy.  When  the  organic  matter  is  completely  broken  down  (oxidized), 
carbon  dioxide  (CO2),  water,  and  minerals  are  left.  The  inorganic  (avail¬ 
able)  form  of  N  resulting  from  decomposition  is  NH'^ . 

The  rate  of  decomposition  and  release  of  available  N  depends  on  the 
nature  of  the  N  compounds,  which  are  greatly  influenced  by  the  type  of 
treatment  or  stabilization  process  the  biosolids  receive,  the  duration  of  the 
process,  and  the  type  of  organic  matter  in  the  biosolids.  Decomposition  is 
slower  when  the  biosolids  are  more  stable.  Mineralization  amounts  are 
higher  in  soils  with  higher  amounts  of  organic  matter;  therefore,  taking 
steps  to  maintain  or  increase  soil  OM  can  help  supply  a  relatively  constant 
amount  of  available  N  to  the  soil.  The  amount  of  mineralization  is  also 
dependent  on  the  type  of  organic  matter  present.  Fresh  manure  or  crop 
residue  will  break  down  faster  than  humus  that  is  the  result  of  years  of 
decomposition.  In  addition,  the  ratio  between  total  soil  carbon  and  total 
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soil  N  affects  how  quickly  this  process  occurs,  because  micro-organisms, 
like  plants,  need  N  to  live. 

Often  N  released  from  biosolids  OM  may  be  taken  up  by  soil  microorgan¬ 
isms  and  converted  back  to  organic  forms.  This  process  is  called  “immobi¬ 
lization.”  N  immobilization  refers  to  the  process  where  inorganic  N  ( N0~ 

or  NH^ )  is  biologically  converted  to  organic  N,  and  is  essentially  the 

reverse  of  N  mineralization.  Immobilization  generally  occurs  in  nutrient- 
poor  soils,  in  soils  with  a  lot  of  matter  that  is  low  in  N  and  high  in  carbon 
(such  as  woody  material  or  straw),  or  in  soils  where  organic  amendments 
(such  as  sawdust  or  low-N  compost)  have  been  added.  NH'^  immobiliza¬ 
tion  appeared  to  be  dominant  when  NO~  and  NH'^  were  both  present. 

It  has  long  been  recognized  that  net  N  mineralization,  the  balance  between 
N  mineralization  or  immobilization,  is  tightly  coupled  with  the  C  cycle. 

The  carbon-to-nitrogen  (C:N)  ratio  is  often  used  as  an  indication  of 
whether  mineralization  or  immobilization  will  occur.  The  C:N  ratio  is  the 
total  concentration  of  C  divided  by  the  total  concentration  of  N.  When  sur¬ 
face  soil  layers  have  a  C:N  ratio  greater  than  30:1,  then  immobilization  is 
highly  likely  to  occur.  This  is  because  microorganisms  need  N  to  assimilate 
the  available  C.  When  the  C:N  ratio  is  below  20:1,  N  mineralization  is 
likely  to  occur.  When  the  C:N  ratio  is  between  20:1  and  30:1,  both 
mineralization  and  immobilization  may  occur  but  they  will  generally 
balance. 

Volatilization 


Ammonia  (NH3)  volatilization  refers  to  the  loss  of  NH^  as  a  gas  into  the 

atmosphere.  The  process  is  increased  at  high  pH  because  NH'^  will  more 

easily  convert  to  NH3  at  high  pH.  The  potential  for  loss  of  NH3  gas  to  the 
air  depends  on  several  equilibrium  relationships  in  the  soil.  NH'^  plays  an 

important  role  in  these  relationships.  Most  of  the  NH'^  is  bound  to  soil 

surface  exchange  sites;  the  remainder  is  dissolved  in  soil  solution.  If  some 
NH'^  is  removed  from  the  soil  solution  through  plant  uptake  or  other 

means,  NH'^  will  move  from  the  soil  surface  sites  into  the  soil  solution  to 
reestablish  equilibrium. 
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The  amount  of  NH^  that  volatilizes  depends  on  a  number  of  environmental 
factors  and  biosolids  management  techniques.  The  important  ones,  in 
addition  to  the  pH,  include  temperature  and  wind  speed,  and  whether  the 
biosolids  are  incorporated,  injected,  or  surface  applied.  As  much  as  lOO 
percent  of  the  initial  NH3  and  NH'^  in  biosolids  that  are  surface  applied  to 

agricultural  soils  can  be  lost  to  volatilization.  Even  when  the  biosolids  have 
been  mixed  in  with  the  soil,  the  majority  of  the  NH3  may  volatilize  if  the 
soil  has  a  high  pH  (over  8.0).  Little  volatilization  may  occur  in  sites  such 
as  forest  environments  because  of  the  low  pH  of  the  forest  floor,  the  low 
wind  speed  in  the  forest  stands,  and  the  low  amounts  of  radiation  reaching 
the  forest  floor. 

Nitrification 

Soil  NH'^  can  quickly  (hours  to  weeks)  be  converted  into  NO^  and  then 

into  NO~ .  NO^  is  an  intermediate  product  in  many  N  transformations. 

This  process,  known  as  nitrification,  only  occurs  in  the  presence  of  oxygen, 
so  generally  it  will  be  slow  or  non-existent  in  waterlogged,  anaerobic  soils. 
During  the  process,  the  microorganism  Nitrosomonas  oxidizes  NH'^  to 

NO^ ,  and  the  microorganism  Nitrobacter  oxidizes  the  NO^  to  NO~ . 

Nitrification  rates  are  highest  when  soils  are  warm  and  moist  and  the  pH 
is  neither  strongly  acidic  nor  alkaline.  Over-application  of  inorganic  N  in 
arid  climates  will  promote  buildup  of  excess  NO~  in  the  soil. 

Denitrification 

Denitrification  refers  to  the  process  where  NO~  becomes  nitrogen  gas 

(N2).  It  is  the  opposite  of  nitrification  in  that  oxygen  is  removed  rather 
than  added.  Denitrification  requires  the  absence  of  oxygen,  or  ‘anaerobic’ 
conditions.  Similar  to  nitrification,  microorganisms  are  responsible  for 
denitrification,  and  therefore  it  occurs  faster  in  warm,  moist  soils.  There¬ 
fore,  denitrification  losses  of  N  are  most  significant  when  soils  alternate 
between  aerobic  conditions,  which  allow  NO~  to  accumulate,  and  anaero¬ 
bic  conditions. 

The  amount  of  denitrification  depends  on  the  availability  of  NO~ ,  soil 

saturation,  soil  temperature,  and  availability  of  easily  decomposable 
organic  matter.  When  a  site  within  the  soil  becomes  saturated  with  water. 
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O2  can  no  longer  diffuse  readily  into  that  site.  If  microorganisms  are  active 
within  that  site,  they  will  deplete  the  O2  and  then  begin  using  the  NO~ .  If  a 

soil  is  too  cold  for  microbial  activity  or  if  there  is  too  little  food  (organic 
matter)  available  for  the  microorganisms,  the  O2  will  not  be  depleted  and 
denitrification  will  not  occur. 

Plant  uptake 

Inorganic  forms  of  N  ( NH^  and  NO~ )  are  taken  up  by  plants  rooted  in  the 

soil  or  floating  in  the  water.  Plant  N  uptake  varies  widely,  depending  on 
plant  type,  growing  conditions  (moisture  and  temperature),  and  manage¬ 
ment  practices.  In  addition,  NO~  requires  more  energy  by  the  plant  after 

it  is  taken  up  because  the  NO~  must  be  converted  to  NH'^  in  the  plant 

before  it  is  made  into  proteins.  It  is  increasingly  being  shown  that  the 
uptake  of  organic  N  can  represent  a  substantial  fraction  of  total  plant  N 
uptake  in  soil  (Bardgett  et  al.  2003). 

Sorption 

Most  soil  profiles  are  negatively  charged  and  sorption  of  N  in  the  soil 
occurs  through  the  process  of  cation  exchange,  whereby  the  NH'^  is 

weakly  bound  to  soil  particles  by  electrostatic  attraction.  The  sorption  pre¬ 
vents  NH^  from  moving  very  rapidly  through  the  soil.  Although  it  may 

seem  that  the  NH^  would  not  be  available  for  plant  uptake,  NH'^  can 

move  away  from  the  soil  surface  as  NH'^  levels  decrease  in  soil  solution 

due  to  the  process  known  as  diffusion.  Hence,  there  is  an  exchange  of 
NH'^  between  soil  and  soil  solution.  More  NH'^  is  held  by  high  pH  (neu¬ 
tral  to  alkaline)  soils,  and  conversely,  NH'^  moves  more  readily  in  low  pH 
(acidic)  soils. 

The  two  negatively  charged  N  forms  ( NO~  and  NO^ )  will  be  repelled  from 

negative  charges  on  the  clay  surface,  and  are  not  attracted  very  strongly  to 
the  fewer  positive  charges  on  clay  surfaces.  Therefore,  NO~  and  NO^ 

both  have  relatively  high  mobility,  meaning  they  can  move  easily  through 
the  soil  and  do  not  undergo  much  exchange. 
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Surface  runoff,  erosion,  and  leaching 

Nitrogen  may  be  lost  from  the  soil  by  surface  runoff  and  erosion.  Dissolved 
forms  of  N  can  be  transferred  from  surface  water  to  soil  solution  (pore- 
water)  and  vice  versa,  through  the  process  of  diffusion.  Erosion  losses  of  N 
can  be  significant.  NO~  is  highly  mobile  as  discussed  previously  and  sub¬ 
ject  to  leaching  losses  when  both  soil  NO~  content  and  water  movement 
are  high.  Movement  of  NO~  through  soil  is  governed  by  convection,  or 

mass  flow,  with  the  moving  soil  solution  and  by  diffusion  within  the  soil 
solution  (Jury  and  Nielson  1989).  NO~  leaching  from  field  soils  must  be 

carefully  controlled  because  of  the  serious  impact  that  it  can  have  on  the 
groundwater.  Leaching  of  NH'^  is  usually  insignificant. 

Soil  nitrogen  mass  balance  equations 

The  basic  concept  in  the  soil  N  module  is  the  mass  balance  that  was 
employed  to  formulate  mathematical  equations  for  each  of  the  soil  N  state 
variables.  As  classified  in  the  SWAT  model,  the  state  variables  of  the  soil  N 
are:  the  fresh  organic  nitrogen  (orgN/rs),  the  active  organic  nitrogen 
(orgNact),  the  stable  organic  nitrogen  (orgNsta),  the  ammonium  N  (NH^), 

and  the  nitrate  N  ( NO~ ).  Fresh  organic  N  is  associated  with  crop  residue 

and  microbial  biomass  while  the  active  and  stable  organic  N  pools  are 
associated  with  the  soil  humus.  NH'^  is  partitioned  between  the  aqueous 

and  sorbed  phases.  Ion  exchange  may  be  an  important  reaction  for  NH'^ 
in  soils,  sediments,  and  aquifer  systems  (Drever  1982).  However,  NO~  is 
very  soluble  in  water,  so  the  corresponding  sorption  of  NO~  is  rare.  There¬ 
fore,  there  are  in  total  six  soil  N  state  variables.  Figure  2  is  a  diagram  of  the 
soil  N  transformation  module  used  in  NSM. 

N  transformation  processes  in  the  soil  are  simulated  in  NSM  for  the  five 
state  variables:  orgNfrs,  orgNact,  orgNsta,  NH'^  and  NO~.  Exchange  of 

NH'^  between  the  soil  solution  and  adsorbed  phase  is  assumed  to  be 

reversible,  NH^  on  soil  particles  is  not  considered  as  a  state  variable  in 

this  NSM  version.  All  N  transformation  rates  are  estimated  using  first- 
order  kinetic  equations  taking  into  account  the  effect  of  soil  water  content 
and  temperature.  This  section  summarizes  the  mass  balance  equations 
used  to  simulate  the  N  cycle  in  soils. 
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Fresh  organic  nitrogen 

This  pool  contains  all  the  crop  N  returned  to  the  soil  during  the  growth  of 
the  crop.  This  includes  N  in  dead  roots,  root  exudates,  fallen  leaves  and 
other  plant  debris  shed  during  the  growing  season.  Upon  decomposition, 
the  N  in  the  crop  residues  will  enter  the  humus  and  biomass  pools.  The 
mass  balance  equation  for  the  soil  fresh  organic  N  concentration  over  time 
is  described  as: 


d(Az-orgNfJ 

dt 


^^min\imb  ^^dec  ^^frs,e  ^^frs,s 


(1) 


where 

orgNfrs  =  concentration  of  soil  layer  fresh  organic  N  pool  [M/L3] 

Az  =  depth  of  the  soil  layer  [L] 

ONmin\imb  =  net  mineralization/immobilization  rate  of  soil  layer  fresh 
organic  N  pool  [M/LVT] 

ONdec  =  decomposition  rate  of  soil  layer  fresh  organic  N  pool  [M/L^/T] 
ONfrs.e  =  net  surface  erosion/deposition  rate  of  soil  layer  fresh  organic 
N  pool  [M/LVT] 

ONfrs.s  =  external  sources  added  to  the  soil  layer  fresh  organic  N  pool 
[M/L2/T]. 

Mineralized  orgNfrsh  is  added  to  the  pool  while  decomposed  orgN/rsh 

is  added  to  the  humus  active  organic  pool,  which  is  further  partitioned  into 
orgNsta  and  orgNact  using  partitioning  ratios  determined  for  soil  organic 
N.  Both  transformation  equations  are  adopted  from  the  SWAT  model.  The 
N  mineralization  rates  in  SWAT  are  net  mineralization  rates,  which 
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incorporate  immobilization  into  the  equations.  N  mineralization  from  the 
residue  fresh  organic  N  is  estimated  as: 


ON 


mm\imb 


Az 


(2) 


where 

fra  =  fraction  of  soil  layer  fresh  organic  N  mineralized  to  the  NH^ 
pool  (0.8) 

Sntr  =  residue  decay  rate  constant  [i/T] 

The  decomposition  equation  is  adopted  from  the  SWAT  model.  Decompo¬ 
sition  from  the  residue  fresh  organic  N  pool  to  the  active  organic  pool  is 
estimated  as: 


ON 


dec 


1 

Az 


(3) 


where 

frh  =  fraction  of  soil  layer  fresh  organic  N  decomposed  to  the  humus 
pool  (0.2) 

Decomposition  and  mineralization  are  controlled  by  a  decay  rate  constant. 
The  decay  rate  constant  is  a  function  of  the  C:N  ratio,  the  C:P  ratio,  and 
the  composition  of  the  crop  residue,  temperature,  and  soil  water. 

^ntr  =  Prsd  ’  Y ntr  (/tmp  ’  Ysw  (4) 


where 

Prsd  =  mineralization  rate  coefficient  of  the  soil  residue  fresh  organic 
N  pool  (0.05  day-1)  [i/T] 

Yntr  =  residue  composition  correction  factor  for  the  soil  layer 
Ytmp  =  temperature  correction  factor  for  the  soil  layer  (>  o.i) 

Ysw  =  water  moisture  correction  factor  for  the  soil  layer  (>  0.05). 

The  residue  composition  correction  factor  is  determined  by: 
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7ntr=min 


exp 


exp 


-0.693 


-0.693 


(^C:N 

25 

(^c:p -200) 


200 


1.0 


(5) 


where 

sc;N  =  C:N  ratio  of  the  residue  in  the  soil  layer 
sc.p  =  C:P  ratio  ofthe  residue  in  the  soil  layer 


In  future  NSM  versions,  C:N  and  C:P  ratios  will  be  calculated  in  the  carbon 
cycle  module  and  passes  to  nutrient  cycling  residue  composition  factor 
equations. 


The  temperature  correction  factor  is  determined  by: 

/tmo  =0-9 - - 7  +  0-1 

r3„,,  +  exp(9.93-0.312-r3,,) 


(6) 


where 

Tsoii  =  temperature  of  soil  layer  [°C] 

The  water  moisture  correction  factor  is  determined  by: 

9 

rsw=^  (7) 


where 


9  =  soil  water  content  (mm  H2O) 

9fc  =  soil  water  content  at  field  capacity  (mm  H2O) 
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Active  organic  nitrogen 

N  enters  this  pool  from  crop  debris,  the  soil  microbial  biomass,  and  from 
any  applied  organic  manure.  The  mass  balance  equation  for  the  soil  active 
organic  N  concentration  over  time  is  described  as: 

c/(Az.orgN.„)  ^ ^ 


where 

orgN^^^  =  concentration  of  soil  layer  active  organic  N  pool  [M/L3] 

ONtrn  =  rate  transferred  between  the  active  and  stable  organic  pool 
[M/L2/T] 

ONmin  =  mineralization  rate  of  soil  layer  active  organic  N  pool  [M/L^/T] 
ONact,e  =  net  surface  erosion/deposition  rate  of  soil  layer  active  organic 
N  pool  [M/LVT] 

ONact,s  =  external  sources  added  to  the  soil  layer  active  organic  N  pool 
[M/L2/T] 


The  organic  N  associated  with  humus  is  partitioned  into  two  pools  (active 
and  stable)  to  account  for  the  variation  in  availability  of  humic  substances 
to  mineralization.  The  active  and  stable  organic  N  pools  are  in  dynamic 
equilibrium  and  are  described  using  a  procedure  adapted  from  the  SWAT 
model.  Organic  N  flux  between  the  active  and  stable  pools  is  estimated  as: 


=  Ptrn  — 

Az 


f  1  ^ 

orgN,,, 

1 

1  fract  ; 

-orgN,,, 

(9) 


where 

yftrn  =  ratc  coustaut  (o. 00001  day  ^  [i/T] 

fract  =  fraction  of  soil  layer  humic  organic  N  in  the  active  organic  N 
pool  (0.02) 

orgNsta  =  concentration  of  soil  layer  stable  organic  N  pool  [M/L3] 

In  SWAT,  N  mineralized  from  the  humus  active  organic  pool  is  added  to 
the  NH^  pool  in  the  layer.  If  the  soil  water  content  exceeds  field  capacity, 

or  soil  temperature  is  less  than  o  °C,  mineralization  from  the  humus  active 
organic  pool  (ammonifi cation)  does  not  occur.  It  is  estimated  as: 
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{/t^p  -/sw)  orgN^ 


(10) 


where 

=  mineralization  rate  coefficient  of  the  soil  humus  active  organic 
N  (0.0003  day-1) 

Stable  organic  nitrogen 

Most  soil  N  retained  in  this  pool  can  accumulate  rapidly  but  is  not  readily 
accessible  to  microbial  mineralization  (Kaye  et  al.  2003).  The  mass  bal¬ 
ance  equation  for  the  soil  stable  organic  N  concentration  over  time  is 
described  as: 


djAz-orgN^J 

dt 


-  O^c/ec  +  ^^trn  “O^sta.e  +  O^sta.s 


(11) 


where 

ONsta,e  =  net  surface  erosion/deposition  rate  of  soil  layer  stable  organic 
N  pool  [M/LVT] 

ONsta,s  =  external  sources  added  to  the  soil  layer  stable  organic  N  pool 
[M/L2/T] 

Ammonium  nitrogen 

Nitrogen  enters  the  soil  NH'^  pool  by  the  mineralization  of  soil  organic  N 

(organic  manure,  humus,  soil  microbial  biomass,  and  crop  debris),  by  the 
hydrolysis  of  urea,  from  any  applied  organic  manure,  and  from  NH^ 

fertilizers.  NH^  is  removed  via  volatilization,  nitrification,  and  by 

immobilization.  NH'^  is  available  to  plants  and  is  not  leached  to  a  great 

extent.  Since  NH^  is  a  positively  charged  ion  (cation),  it  may  be  attracted 

to  and  held  by  the  negatively  charged  soil  clay.  The  mass  balance  equation 
for  the  soil  NH'^  concentration  over  time  is  described  as: 

-  NH„„„  -  .  +  NH,  ( 12) 
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where 

NH^  =  concentration  of  soil  layer  NH^  pool  [M/L3] 

NHmin  =  total  mineralization  processes  rate  of  soil  layer  organic  N  pools 
[M/L2/T] 

NHnit\voi  =  net  nitrification/volatilization  processes  rate  in  the  soil  layer 
[M/L2/T] 

NHup  =  plant  uptake  rate  of  soil  layer  NH'^  pool  [M/L^/T] 

RNH4,e  =  mass  transfer  rate  of  NH'^  between  the  upper  soil  layer  and 
surface  runoff  [M/L^/T] 

NHs  =  external  sources  added  to  the  soil  layer  NH^  pool  [M/L^/T] 

Based  on  Equations  2  and  10,  the  total  mineralization  rate  to  the  NH'^ 
pool  is  calculated  as: 

1  1  12 

Nitrification  and  volatilization  equations  are  adopted  from  the  SWAT 
model.  The  total  rate  of  lost  to  nitrification  and  volatilization  is  cal¬ 
culated  using  a  first-order  kinetic  rate  equation  (Reddy  et  al.  1979): 

^  ^^4"  [1  -  exp  (-77™t  -  %o> )]  (14) 


where 


rjnit  =  nitrification  regulator 
rjvoi  =  volatilization  regulator 

To  partition  NHnit\voi  between  nitrification  and  volatilization,  the  expres¬ 
sion  by  which  is  multiplied  in  Equation  14,  is  solved  using  each  regu¬ 
lator  individually  to  obtain  a  fraction  of  each  process.  The  rate  of  N 
removed  from  the  pool  by  nitrification  is  then  calculated  as: 


= 


{Kit  +  Ki) 


NH. 


nit\vol 


(15) 
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and  the  rate  of  N  removed  from  the  NHl  pool  by  volatilization  becomes: 


fr.„ 


[Kt  +  K,) 


NH. 


nit\vol 


(16) 


where 


NHnit  =  nitrification  processes  rate  in  the  soil  layer  [M/L^/T] 

NHvoi  =  volatilization  processes  rate  in  the  soil  layer  [M/L^/T] 
fruit  =  i-exp{-rinit)  (estimated  fraction  of  N  lost  by  nitrification) 
frvoi  =  i-exp{-T]voi)  (estimated  fraction  of  N  lost  by  volatilization) 

The  impact  of  environmental  factors  on  nitrification  and  ammonia  volatili¬ 
zation  is  defined  by  the  nitrification  regulator  and  volatilization  regulator. 
The  nitrification  regulator  is  determined  by: 


7n/t  =7tmp-7sw  (17) 

and  the  volatilization  regulator  is  determined  by: 

rival  ='ntmp-r!voiz  (18) 


where 

Jjtmp  =  nitrification/volatilization  temperature  correction  factor 
rjsiv  =  nitrification  soil  water  correction  factor 
rjvoiz  =  volatilization  depth  correction  factor 

The  nitrification/volatilization  temperature  correction  factor  is  estimated 
as: 


77  =0  41^^^°'' _ ^  if  T  >5 

Vtmp  II  ‘soil^^ 


(19) 


The  nitrification  soil  water  correction  factor  is  estimated  as: 


7sw  = 


O.25{0,,-0^,) 


if  0-0^,<O.25(0,a-0^p) 


(20) 
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73.  =1.0  if  0-0„,>O.25{0,,-0^,) 


where 


9wp  =  soil  water  content  at  wilting  point  (mm  H2O) 
The  volatilization  depth  correction  factor  is  estimated  as: 


rivoiz  =  1 


z 

z +  exp(4.706 -0.305 -z) 


(21) 


(22) 


where 


z  =  depth  from  the  soil  surface  to  the  middle  of  the  layer  [L] 

Plant  uptake  is  a  very  complicated  process,  involving  crop,  environment, 
and  cultivation  factors.  The  uptake  rate  is  known  to  depend  both  on  a 
plant’s  requirements  and  the  availability  of  N  in  the  soil  to  meet  those 
requirements.  The  crop  N  uptake  is  limited  by  the  inorganic  N  available  in 
the  soil  layer.  Adopted  from  the  NITS  model  (Birkinshaw  and  Ewen 
2000),  the  rate  of  plant  uptakes  from  is  estimated  as: 


/V/-/^p  =  mini 


/VH 


nh: 


up 


/VH;+/V03 


(23) 


where 

=  potential  plant  uptake  rate  of  soil  layer  [M/L^/T] 
kup  =  maximum  plant  uptake  rate  of  soil  layer  [T-i] 

Nitrate  nitrogen 

N  enters  the  soil  NO~  pool  from  NO~  fertilizers,  inputs  from  the  atmos¬ 
phere,  and  from  the  nitrification  of  NH'^ .  NO~  is  one  of  the  principal 
forms  of  N  used  by  plants.  Unlike  NH'^ ,  NO~  exists  only  in  the  dissolved 

phase  and  is  very  soluble,  being  subject  to  only  dissolved  transport,  and 
moves  through  the  soil  profile  with  percolating  water.  In  modeling  NO~ , 
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the  partition  coefficient  equals  zero.  The  mass  balance  equation  for  the  soil 
NO~  concentration  over  time  is  described  as: 


djAz-NQ-^) 

dt 


“  ^^nit  ^^dnit  ^N03,f  ^N03,e 


(24) 


where 

NO~  =  concentration  of  soil  layer  NO~  pool  [M/L3] 

NOdnit  =  denitrification  processes  rate  in  the  soil  layer  [M/LVT] 
NOup  =  plant  uptake  rate  of  soil  layer  NO~  pool  [M/L^/T] 

RN03,e  =  mass  transfer  rate  of  NO~  between  the  upper  soil  layer  and 
surface  runoff  [M/L^/T] 

Rno3/  =  infiltration  rate  of  soil  layer  NO~  pool  [M/LVT] 

NOs  =  external  sources  added  to  the  soil  layer  NO~  pool  [M/L^/T]. 


The  iVOg  form  of  N  is  a  major  concern  in  groundwater  pollution.  Leaching 
of  NO~  from  the  soil  layers  is  modeled  as: 


^N03,f  ~  ^  ■  ^^3 


(25) 


Denitrification  of  NO~  to  N2  is  assumed  to  occur  when  soil  water  exceeds 
field  capacity.  Based  on  the  SWAT  model,  the  denitrification  rate  of  NO~ 
is  estimated  as: 


NO,,,  =N0; 


1 -exp  (-1.4- 7, „p-orgC)]  if  73w^0.95 


(26) 


/VO,,,  =0.0  if  73w<0.95 


(27) 


Adopted  from  the  NITS  model  (Birkinshaw  and  Ewen  2000),  the  rate  of 
plant  uptakes  from  NO~  is  estimated  as: 


N0,p  =min<l 


NO, 


NO, 


up 


/VO3 +/vh; 


Kno, 


(28) 
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where 

NOlp  =  potential  plant  uptake  rate  of  soil  layer  NO~  [M/L^/T] 
kup  =  maximum  plant  uptake  rate  of  soil  layer  NO~  [i/T] 

In  future  NSM  versions,  plant  uptake  of  N  will  be  calculated  in  the  plant 
dynamic  module  and  passed  to  the  inorganic  N  mass  balance  equations 
described  above. 

Initialization  of  nitrogen  levels 

In  addition  to  the  appropriate  transport  fields  given  by  hydrologic  and 
sediment  transport  simulations,  the  user  must  define  the  initial  N  levels 
for  the  state  variables  within  the  watershed.  The  following  is  provided  as 
an  empirical  reference  if  measured  data  are  not  available. 

Nitrogen  in  rainfaii 

Lightning  discharge  converts  atmospheric  N2  to  nitric  acid,  which  can  then 
be  transferred  to  the  soil  with  precipitation.  It  has  been  estimated  that 
rainfall  adds  about  10  lb  of  N  to  the  soil  per  acre  per  year.  More  N  will  be 
added  to  the  soil  with  rainfall  in  areas  with  a  high  amount  of  lightning 
activity  than  in  areas  with  little  lightning.  In  SWAT,  the  N  in  rainfall  is 
added  to  the  nitrate  pool  in  the  top  10  mm  of  soil.  The  concentration  of 
nitrate  added  to  the  soil  in  rainfall  is  calculated  as: 

(29) 

h 


where 

NOrain  =  coucentratiou  of  NO~  added  by  rainfall  [M/L3/T] 

Rno3  =  concentration  of  N  in  the  rain  [M/L3] 

R  =  precipitation  intensity  (mm  H2O)  [L/T] 
h  =  precipitation  depth  [L] 

Nitrogen  in  soiis 

The  user  will  be  able  to  either  define  the  initial  N  levels  or  use  model 
defaults  to  initialize  the  levels  of  N  in  the  different  pools.  If  the  user  allows 
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the  model  to  initialize  the  N  levels,  then  the  method  as  described  below 
will  be  followed. 


Inorganic  forms  of  N  are  usually  added  to  the  soil  by  precipitation  (i.e., 
rain  or  snowfall),  or  as  fertilizers.  Microorganisms  in  the  soil  convert 
organic  forms  of  N  into  inorganic  forms  that  are  then  usable  by  plants.  The 
ammonium  N  pool  for  the  soil  NH'^  is  initialized  to  o  ppm.  Initial  nitrate 

levels  in  the  soil  are  varied  by  depth  using  the  following  relationship: 


NO 


3,z 


=  7 • exp 


f— 1 

Uoooj 


(30) 


where 


z  =  depth  from  the  soil  surface  [L] 

NO~^  =  concentration  of  NO~  in  the  soil  at  depth  z  [M/M  or  ppm],  the 
NO~  concentration  for  a  layer  is  calculated  by  solving  Equa¬ 
tion  30  for  the  soil  horizon’s  lower  boundary  depth. 

Nitrogen  enters  the  soil  in  organic  forms  such  as  plant  roots,  leaves,  and 
other  plant  materials,  in  addition  to  dead  animals,  insects,  and  microor¬ 
ganisms,  manure,  compost,  and  sewage  sludge.  Organic  N  levels  will  be 
assigned  assuming  that  the  initial  C:N  ratio  for  humic  materials  is  14:1. 
The  concentration  of  humic  organic  N,  in  a  soil  layer,  is  calculated  as: 


orgN 


hum 


=  10 


orgC 

14 


(31) 


where 

orgC  =  concentration  of  organic  carbon  in  the  soil  layer  (%) 
orgNhum  =  concentration  of  soil  humic  organic  N  [M/M  or  ppm] 

The  humic  organic  N  is  partitioned  between  the  active  and  stable  pools 
using  the  following  equations: 


orgN,,t^orgN,^^-fr^^,  (32) 

orgN^,^  =  orgN,^^{l-fr^^,) 


(33) 
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In  SWAT,  the  fraction  of  humic  N  (fract)  in  the  active  N  pool  is  set  to  0.02. 
N  in  the  fresh  organic  pool  will  be  set  to  zero  in  all  layers  except  the  top 
10  mm  of  soil.  In  the  top  10  mm,  the  fresh  organic  N  pool  will  be  set  to 
0.0015  of  the  initial  concentration  of  residue  on  the  soil  surface. 


orgNf^^  =  0.0015  •  rsd 


(34) 


where 

rsd  =  concentration  of  residue  in  the  soil  layer  [M/L3] 

Soil  phosphorus  module 

Phosphorus  cycle 

Phosphorus  can  exist  in  the  soil  as  phosphate  {HPO~^  or  H^PO~),  par¬ 
ticulate  P,  organic  P,  or  in  P  minerals.  Phosphate  is  the  only  form  that 
plants  can  take  up,  yet  in  most  agricultural  soils  there  is  less  than  1  mg/L 
(ppm)  of  phosphate  in  solution,  which  represents  much  less  than  1  percent 
of  the  total  soil  P.  Organic  P,  which  is  P  bound  in  organic  matter,  has  been 
found  to  represent  between  about  25  percent  and  65  percent  of  total  P  in 
surface  soils,  with  mineral  P  (such  as  calcium  phosphate  minerals)  and 
sorbed  P  representing  the  remainder  (Brady  1984). 

Phosphorus  is  generally  much  less  mobile  than  N,  being  strongly  adsorbed 
to  soil  particles  as  well  as  organic  matter.  Figure  3  shows  the  major  com¬ 
ponents  of  the  watershed  model  P  cycle.  Soil  P  is  divided  into  an  organic 
and  mineral  component  that  can  receive  inputs  via  inorganic  fertilizers, 
organic  manure,  waste,  and  sludge.  Organic  P  in  the  soil  is  divided  into 
three  pools:  the  fresh  organic  P  pool  {orgPfrs),  the  active  organic  P  pool 
(orgPact)  and  the  stable  organic  P  pool  (orgPsta).  The  orgPfrs  represents  the 
organic  matter  that  can  be  easily  mineralized  (e.g.  manures,  decayed 
plants,  and  microbial  biomass).  The  orgPact  constitutes  materials  with  a 
slower  rate  of  mineralization,  whereas  orgPsta  is  composed  of  P  in  stable 
organic  matter  (i.e.,  humus).  Organic  pools  are  differentiated  on  the  basis 
of  C:P  ratios.  Inorganic  P  is  divided  into  labile  (soluble)  inorganic  P  (Pso/), 
active  inorganic  P  (minPact),  and  stable  inorganic  P  iminPsta)  pools.  Psoi  is 
in  rapid  equilibrium  with  minPact,  which  in  return  is  in  slow  equilibrium 
with  minPsta-  When  inorganic  fertilizer  P  is  added,  it  rapidly  equilibrates 
between  Psoi  and  minPact.  The  slow  reaction  between  minPact  and  minPsta 
then  follows. 
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P  transformation  processes  in  NSM  are  simulated  for  organic  and  inor¬ 
ganic  P  in  the  soil.  The  processes  that  control  the  concentration  of  avail¬ 
able  P  in  the  soil  are:  mineralization,  immobilization,  precipitation,  disso¬ 
lution,  adsorption,  desorption,  plant  uptake,  surface  runoff,  and  erosion. 
Each  of  these  processes  is  described  below. 

Soil  phosphorus  transformations 

The  P  transformations  in  soils  involve  complex  mineralogical,  chemical, 
and  biological  processes.  P  transformations  in  the  soil  include  decomposi¬ 
tion  and  mineralization  of  organic  P,  immobilization  of  labile  P,  and  sorp¬ 
tion  of  labile  P  to/from  soil  particles.  P  dissolved  in  soil  water  and  P 
adsorbed  to  soil  particles  can  be  transported  to  surface  water.  Microorgan¬ 
isms  deeply  affect  soil  P  transformation  through  mineralization- 
immobilization  processes.  Knowledge  of  P  transformations  is  essential  to 
understanding  P  behavior  in  soils. 
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Mineralization  and  immobilization 

In  general,  P  mineralization  and  immobilization  are  similar  to  those  of  N 
in  that  both  reactions  or  processes  occur  simultaneously  in  soils  and  can 
be  depicted  as  follows: 


Mineralization  converts  organic  P  to  inorganic,  plant-available  P.  Immobi¬ 
lization  is  the  reverse  process  of  mineralization  where  inorganic  P  is  con¬ 
verted  back  to  plant  unavailable  organic  P.  The  C:P  ratio  of  the  decompos¬ 
ing  residues  regulates  the  predominance  of  P  mineralization  over  immobi¬ 
lization,  just  as  the  C:N  ratio  regulates  N  mineralization  over  immobiliza¬ 
tion.  The  following  guidelines  have  been  suggested:  mineralization  occurs 
most  readily  when  the  C:P  ratio  is  less  than  200:1,  and  immobilization 
occurs  when  that  ratio  is  greater  than  300:1  (Havlin  et  al.  1999).  Minerali¬ 
zation  and  immobilization  of  P  are  affected  by  temperature,  moisture, 
aeration,  and  pH  in  similar  ways  as  N  mineralization  and  immobilization, 
because  they  involve  the  same  microbial  processes. 

Adsorption  and  desorption 

Adsorption  refers  to  the  binding  of  P  to  soil  particles.  Because  phosphate 
has  a  negative  charge,  it  is  attracted  to,  and  binds  strongly  to,  positively 
charged  minerals,  such  as  aluminum  (Al)  and  iron  (Fe)  hydroxides  and 
oxides.  Like  other  soil  particles,  these  minerals  become  more  positively 
charged  at  lower  pH;  therefore,  more  phosphate  is  sorbed  at  a  lower  pH. 
Finer  textured  soils  generally  can  adsorb  more  P  because  they  have  more 
surface  area. 


Adsorption  is  decreased,  and  hence  available  P  levels  are  increased,  when 
the  soil  solution  contains  high  levels  of  other  anions  such  as  bicarbonate, 
carbonate,  silicate,  sulfate,  or  molybdate  that  compete  for  sorption  ‘sites.’ 
In  addition,  dissolved  organic  compounds  associated  with  organic  matter 
can  increase  P  availability  by  competing  for  phosphate  sorption  sites  or 
coating  Fe/Al  oxides.  Adsorption  of  P  generally  increases  with  increased 
temperatures  in  P  fertilized  soils,  but  there  is  no  correlation  with  tempera¬ 
ture  in  soils  that  have  not  had  P  added  (Havlin  et  al.  1999).  Desorption  is 
the  opposite  of  adsorption.  P  desorption  generally  increases  as  solution  P 
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decreases,  or  under  flooded  conditions  due  to  dissolution  of  Fe  hydroxides 
and  oxides,  that  release  adsorbed  P. 

Precipitation  and  dissoiution 

Precipitation  is  the  process  through  which  soluble  P  is  converted  to  part  of 
mineral  P.  Available  P  concentrations  are  largely  controlled  by  the  solubil¬ 
ity  of  P  minerals  that  are  dominated  by  calcium  phosphates  (Ca-P)  in  neu¬ 
tral  to  high  pH  soils,  and  by  Al  and  Fe  phosphates  (Al-P  and  Fe-P)  at  pH 
levels  below  about  6.5.  There  are  numerous  forms  of  calcium  phosphates 
in  soil,  ranging  from  the  very  soluble  monocalcium  phosphate  (MCP)  to 
the  very  insoluble  fluorapatite.  After  fertilizing  with  P  in  a  neutral  or  high 
pH  soil,  MCP  will  form  first,  followed  by  the  other  calcium  phosphates  in 
order  from  high  to  low  solubility.  The  time  for  each  mineral  to  form  is 
highly  dependent  on  temperature.  If  a  soil  has  a  mixture  of  the  various  Ca- 
P  solids,  the  more  soluble  forms  will  dissolve  more  readily  as  phosphate 
levels  in  solution  decrease  during  the  growing  season.  Decreasing  pH 
would  take  large  quantities  of  acid,  or  acid-producing  substances  such  as 
elemental  S,  Fe+3,  or  manure.  As  may  be  expected,  soils  with  higher  levels 
of  calcium  carbonate  (lime)  will  tie  up  more  P  due  to  precipitation  of  Ca 
phosphates,  and  thus  lower  the  soil  test  P. 

Al  phosphates  and  Fe  phosphates  are  the  predominant  P  minerals  in  soils 
with  pH  levels  below  about  6.5  (Havlin  et  al.  1999).  The  solubility  of  these 
minerals  decreases  at  lower  pH,  directly  opposite  of  the  solubility  for  cal¬ 
cium  phosphates.  Therefore,  P  is  most  available  around  pH  6.5,  because  at 
lower  pH  levels,  P  retention  is  high  due  to  Al-P  and  Fe-P  precipitation,  and 
at  higher  pH  levels,  Ca-P  minerals  precipitate.  Precipitation  and  dissolu¬ 
tion  are  not  depicted  within  the  current  NSM  soil  P  module. 

Piant  uptake 

Despite  low  concentrations  of  phosphate  in  soil  solution,  plants  can  take 
up  substantial  amounts  of  P  due  to  P  desorption  and  dissolution,  followed 
by  P  diffusion  to  the  plant  root  (Foth  and  Ellis  1997).  By  taking  up  large 
amounts  of  P,  a  strong  ‘diffusion  gradient’  is  created,  which  moves  P 
toward  the  plant  root  at  much  higher  rates  than  water  is  moving  via  tran¬ 
spiration.  Plant  dynamic  modules  are  currently  under  development  and 
will  be  available  in  future  releases  of  NSM. 
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Erosion  and  surface  runoff 

Erosion  represents  a  loss  of  P  from  soils  and  can  occur  from  water  or  wind 
because  P  is  generally  bound  so  tightly  to  the  soil.  Dissolved  P  in  runoff 
can  represent  another  loss  of  P.  However,  the  concentration  of  dissolved  P 
in  runoff  is  generally  low  due  to  the  high  amount  of  sorption  and  precipita¬ 
tion  of  P  minerals.  One  exception  to  this  general  rule  is  animal  manure 
stockpiles  or  manure  application  sites,  where  P  is  concentrated  and  sorp¬ 
tion  sites  in  the  manure  and  surrounding  soil  may  be  approaching  satura¬ 
tion.  It  is  important  to  realize  that  the  general  perception  that  P  binds 
strongly  to  soil  is  true  only  to  a  certain  soil  P  level  or  application  amount. 

A  large  amount  of  research  has  correlated  soil  test  P  levels  with  dissolved  P 
in  soil  solution,  drainage  water,  or  runoff.  One  of  these  studies  found  that 
P  in  solution  began  to  significantly  increase  only  when  Olsen  P  levels  were 
above  6o  ppm  (Heckrath  et  al.  1995). 

Soil  phosphorus  mass  balance  equations 

Simulation  of  soil  P  dynamics  is  taken  from  the  SWAT  model  in  which  six 
pools  (state  variables)  of  P  are  identified.  They  are:  the  fresh  organic  P 
(orgPfrs),  the  active  organic  P  {orgPact),  the  stable  organic  P  (orgPsta),  the 
active  inorganic  P  {minPact),  the  stable  or  insoluble  inorganic  P  {minPsta), 
and  the  labile  (soluble)  P  (PsoO-  Figure  4  is  a  diagram  of  the  soil  P  transfor¬ 
mation  module  used  in  NSM. 


Inorganic  firtilizer 
or  air  deposition 


Erosion  Deposition 


Mass  transfer 


Plant  uptake 


Erosion 


Deposition 


Organic  source: 
manue,  residue 


Soils 


Adsorption 


Leaching 


Inorganic  Pool 


Mineralization 


Figure  4.  Soil  phosphorus  transformation  module. 


P  transformation  processes  in  the  soil  are  simulated  by  NSM  for  the  six 
pool  state  variables:  orgPfrs,  orgPact,  orgPsta,  Psoi,  minPact,  and  minPsta-  All 
P  transformation  rates  are  estimated  using  first-order  kinetic  equations 
taking  into  account  the  effect  of  soil  water  content  and  temperature.  This 
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section  summarizes  the  mass  balance  equations  used  to  simulate  the  P 
cycle  in  soils. 

Fresh  organic  phosphorus 

This  pool  consists  of  P  in  undecomposed  plant  residues,  livestock  excre¬ 
tion,  manure  and  microbes.  The  mass  balance  equation  for  the  soil  fresh 
organic  P  concentration  over  time  is  described  as: 

d(Az  ■  orgPfJ  QP  _Qp  _Qp  Qp  pg. 

^‘min\imb  ^‘dec  ^‘frs,e^^‘frs,s 

where 

orgPfrs  =  concentration  of  soil  layer  fresh  organic  P  pool  [M/L3] 

OPdec  =  decomposition  rate  of  soil  layer  fresh  organic  P  pool  [M/LVT] 
OPminiimb  =  net  mineralization/immobilizatiou  rate  of  soil  layer  fresh 
organic  P  pool  [M/LVT] 

OPfrs,e  =  net  surface  erosion/deposition  rate  of  soil  fresh  organic  P  pool 
[M/L2/T] 

OPfrs,s  =  external  sources  added  to  the  soil  layer  fresh  organic  P  pool 
[M/L2/T] 

Mineralized  orgPfrsh  is  added  to  soluble  P  while  decomposed  orgPfrsh  is 
added  to  humus  organic  pool  (orgPhum).  The  orgPhum  is  further  partitioned 
into  orgPsta  and  orgPact  using  partitioning  ratios  determined  for  soil 
organic  P.  Both  transformation  equations  are  adopted  from  the  SWAT 
model.  The  P  mineralization  rates  are  net  mineralization  rates  that  incor¬ 
porate  immobilization  into  the  equations.  P  mineralized  from  the  fresh 
organic  pool  is  added  to  the  soluble  P  pool  and  is  estimated  as: 

Sntr  •  orgPf,,  (37) 


Phosphorus  decomposed  from  the  fresh  organic  pool  is  added  to  the 
humus  organic  pool.  It  is  estimated  as: 


(38) 
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/Act/Ve  organic  phosphorus 

This  pool  is  P  in  the  solid  phase,  which  is  relatively  easily  transformed  to 
the  soil  solution.  The  mass  balance  equation  for  the  soil  humic  active 
organic  P  concentration  over  time  is  described  as: 

d(AZ.OfgP„)  ^  ^  jggj 


where 

orgPact  =  concentration  of  soil  layer  active  organic  P  pool  [M/L3] 

OPmin  =  mineralization  rate  of  soil  humic  active  organic  P  pool 
[M/L2/T] 

OPtrn  =  rate  transferred  between  the  active  and  stable  organic  P  pools 
[M/L2/T] 

OPact,e  =  net  surface  erosion/deposition  rate  of  soil  humic  active 
organic  P  pool  [M/LVT] 

OPact.s  =  external  sources  added  to  the  soil  layer  active  organic  P  pool 
[M/L2/T] 

P  mineralized  from  the  humus  active  organic  pool  is  added  to  the  soluble  P 
pool  in  the  layer.  It  is  estimated  as: 

^  •  ^’nin  (rtmp  '  /sw  OrgP^^^  (40) 


Stable  organic  phosphorus 

The  soil  stable  organic  P  is  relatively  inert  to  transformation.  The  mass 
balance  equation  for  the  soil  stable  organic  P  concentration  over  time  is 
described  as: 


djAz-orgP^^ 

dt 


=  OPdec+ OPtrn -OPsta.e+OPsta.s 


(41) 


where 

orgPsta  =  concentration  of  soil  layer  stable  organic  P  pool  [M/L3] 
OPsta,e  =  net  surface  erosion/deposition  rate  of  soil  humic  stable 
organic  P  pool  [M/LVT] 


ER  DC/EL  TR-06-12 


32 


OPsta,s  =  external  sources  added  to  the  soil  layer  stable  organic  P  pool 
[M/L2/T]. 

Soluble  phosphorus 

The  soil  soluble  P  pool  is  very  small  and  will  usually  contain  only  a  fraction 
of  a  pound  of  P  per  acre.  The  soluble  P  usually  is  in  the  inorganic  form,  but 
small  amounts  of  organic  P  may  exist  as  well.  Of  total  soluble  P,  organic  P 
is  more  labile  than  inorganic  P  since  it  is  not  easily  bound  to  soil  particles. 
Inorganic  P  is  a  preferred  form  for  plant  growth  although  organic  P  is  also 
uptaken.  The  inorganic  P  in  solution  not  uptaken  by  plants  or  immobilized 
by  microorganisms  can  be  adsorbed  to  mineral  surface  or  precipitated  as 
secondary  P  compounds.  With  low  soluble  P  concentrations,  adsorption 
appears  to  dominate,  while  precipitation  reactions  proceed  when  the  con¬ 
centration  of  P  and  associated  cations  in  the  soil  solution  exceeds  that  of 
the  solubility  product  of  the  mineral.  The  mass  balance  equation  for  the 
soil-soluble  P  concentration  over  time  is  described  as: 

='Pn,i„  -'P»  +'P.  (42) 

where 

Psoi  =  concentration  of  soil  layer  soluble  P  pool  [M/L3] 

IPmin  =  total  mineralization  processes  rate  of  soil  layer  organic  P  pools 
[M/L2/T] 

IPsoi\act  =  net  sorption/desorption  rate  transferred  between  the  soluble  P 
pool  and  active  inorganic  P  pool  [M/L^/T] 

IPup  =  plant  uptake  rate  of  soil  layer  soluble  P  pool  [M/L^/T] 

RDip,e  =  mass  transfer  rate  of  soluble  P  between  the  upper  soil  layer 
and  surface  runoff  [M/L^/T] 

IPs  =  external  sources  added  to  the  soil  layer  soluble  P  pool 
[M/L2/T] 

Based  on  Equations  37  and  40,  the  total  mineralization  rate  to  the  soluble 
P  is  calculated  as: 

IPmin  =  ^  {rtmp  '  /sw  ^  ^  (43) 
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Equilibrium  sorption  between  the  solution  and  active  inorganic  P  pool  is 
governed  by  the  P  availability  index  (PAI).  Sharpley  et  al.  (1984)  have  esti¬ 
mated  values  of  the  PAI  for  different  soils  based  on  regression  analyses. 
PAI  is  the  P  sorption  coefficient  defined  as  the  fraction  of  fertilizer  P 
remaining  in  the  labile  pool  after  the  initial  rapid  phase  of  P  sorption  is 
complete.  The  PAI  is  determined  by: 


PAI  = 


P  -P 

'sol,f  ‘sol,i 


soijrt 


(44) 


where 

Psohf  =  concentration  of  P  in  solution  after  fertilization  and  incubation 
[M/L3] 

Psoij  =  concentration  of  P  in  solution  before  fertilization  [M/L3] 

Psohfrt  =  concentration  of  soluble  P  fertilizer  added  to  the  soil  [M/L3] 


SWAT  simulates  P  sorption  by  assuming  the  soluble  inorganic  P  pool  is  in 
rapid  equilibrium  with  the  active  inorganic  P  pool.  The  transformation  of 
P  between  the  solution  and  active  inorganic  pools  is  determined  by  the  fol¬ 
lowing  equations: 


=  — 


sol\act 


Az 


Pso,-mi^Paot 


PAI 

1-PAI 


If  P^^,>mmP^^^ 


f  PAI 
yl-PAI 


(45) 


IP  = _ B 

so/|act  Htm 


P,o,-minP^,, 


PM 

1-PAI 


If  <  minP^^, 


PAI 

1-PAI 


(46) 


where 

/Jtm  =  first  order  rate  constant  for  optimal  temperature  and  moisture 
conditions  (0.1  day  O 

Uptake  of  P  (mainly  in  the  form  oiH^PO^ )  by  plant  roots  is  assumed  to  fol¬ 
low  Michaelis-Menten  kinetics  (Barber  1980)  and  is  limited  by  soil  water 
availability: 


IP... 


sol 


up 


Mju  +  P,,i 


(47) 
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where 

Umax  =  maximum  rate  of  P  uptake  by  the  given  plant  [M/L^/T] 

/I  =  ratio  of  potential  crop  coefficient  for  a  given  time-step  to  the 
maximum  crop  coefficient  i=Kco/ Kco(max)) 

M  =  Michaelis-Menten’s  half-saturation  constant  of  soluble  P  con¬ 
centration  for  the  given  plant  [M/L3] 
k(p  =  soil  water  stress  factor 

The  soil  water  stress  factor  is  applied  to  limit  the  P  uptake  rate  according 
to  soil  water  availability.  There  is  no  limitation  when  the  soil  water  content 
exceeds  the  field  capacity,  but  P  uptake  decreases  linearly  as  the  soil  water 
content  decreases  from  field  capacity  to  wilting  point  and  no  uptake  occurs 
when  the  soil  water  content  is  below  the  wilting  point.  Thus,  Equation  47 
reflects  both  the  P  demand  and  availability  in  the  root  zone.  In  future  NSM 
versions,  crop  uptake  P  will  be  calculated  in  the  plant  dynamic  module  and 
passed  to  the  inorganic  P  mass  balance  equation. 

Active  inorganic  phosphorus 

Active  inorganic  P,  which  is  assumed  to  be  adsorbed  on  soil  particles, 
represents  more  stable  P  that  is  not  easily  desorbable  but  is  in  equilibrium 
with  soluble  inorganic  P.  As  plants  take  up  phosphate,  the  concentration  of 
phosphate  in  solution  is  decreased  and  some  phosphate  from  the  active 
inorganic  P  pool  is  released.  Because  the  soluble  P  pool  is  very  small,  the 
active  inorganic  P  pool  is  the  main  source  of  available  P  for  crops.  The 
mass  balance  equation  for  the  soil  active  inorganic  P  concentration  over 
time  is  described  as: 


d(Az-minP^J 

dt 


sol\act 


act\sta 


-IP  +  IP 

“  act.e  “  act,s 


(48) 


where 

minPact  =  concentration  of  soil  layer  active  inorganic  P  pool  [M/Ls] 
IPact\sta  =  net  slow  sorption/desorption  transfer  rate  between  the  active 
inorganic  P  pool  and  the  stable  inorganic  P  pool  [M/L^/T] 
IPact,e  =  surface  erosion/deposition  rate  of  soil  active  inorganic  P 
detachment  [M/L^/T] 

IPact,s  =  external  sources  added  to  the  soil  layer  active  inorganic  P  pool 
[M/L2/T] 
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The  active  and  stable  inorganic  P  pools  are  dynamic  and  are  assumed  to 
exist  at  a  slow  equilibrium  with  the  stable  mineral  P  pool.  At  equilibrium, 
the  stable  inorganic  pool  is  four  times  the  size  of  the  active  inorganic  pool. 
The  transformation  of  P  between  the  active  and  stable  inorganic  pools  is 
determined  by  the  following  equations: 

/Pact|sta=^^eqs(4-niinP,,,-minP3,J  If  minP3,3<4-minP,,,  (49) 

± 

^^act\sta  ~  ~7  /^trn  '  Peqs  (4-minP3,,-minP3,3)  If  minP3,3  > 4-minP3,,  (50) 


where 


/3eqs  =  slow  equilibrium  rate  constant  (0.0006  day  i)  [i/T] 


Stable  inorganic  phosphorus 

This  pool  comprises  stable  forms  of  inorganic  P  adsorbed  on  the  soil  in 
equilibrium  with  the  active  pool  of  inorganic  P.  The  inorganic  phosphate 
compounds  in  this  pool  are  more  crystalline  in  their  structure  and  less 
soluble  than  those  compounds  considered  to  be  in  the  active  inorganic  P 
pool.  Some  slow  conversion  between  the  stable  P  pool  and  the  active  P 
pool  does  occur  in  soils.  The  mass  balance  equation  for  the  soil  stable  inor¬ 
ganic  P  concentration  over  time  is  described  as: 


cf(Az-minP3,J 

dt 


=  IP  -IP  +IP 

''  act\sta  “  sta,e  ''  sta,s 


(51) 


where 

minPsta  =  concentration  of  soil  layer  stable  inorganic  P  [M/L3] 

IPsta,e  =  surface  erosion/deposition  rate  of  soil  stable  inorganic  P 
detachment  [M/L^/T] 

IPsta,s  =  external  sources  added  to  the  soil  layer  stable  inorganic  P  pool 
[M/L2/T] 

Initialization  of  soil  phosphorus  levels 

Users  can  define  the  concentration  of  inorganic  and  organic  P  in  all  soil 
layers  at  the  beginning  of  the  simulation.  If  such  information  is  not 
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available,  the  following  default  values  will  be  used  to  initialize  levels  of  P  in 
the  different  pools. 

Initial  soluble  P  pool  for  the  soil  is  set  at  5  mg/kg  to  simulate  unmanaged 
conditions.  A  concentration  of  25  mg/kg  is  considered  representative  of 
cropland.  Organic  P  decreases  quickly  with  soil  depth,  paralleling 
decreases  in  organic  matter.  The  distribution  of  organic  P  with  depth  also 
varies  among  soils.  The  concentrations  of  P  in  the  active  inorganic  pool 
and  the  stable  inorganic  pool  are  initialized,  respectively,  to  (Jones  et  al. 
1984): 


1-PAI 

minP3,=P3, 

(52) 

minP3fa=4-minP,,j 

(53) 

Organic  P  levels  are  assigned  assuming  that  the  N:P  ratio  for  humic  mate¬ 
rials  is  8:1.  The  concentration  of  humic  organic  P  in  a  soil  layer  is 
calculated  as: 


orgp.ur.^  0.125 -orgN,^^  (54) 

P  in  the  fresh  organic  pool  is  set  to  zero  in  all  layers  except  the  top  10  mm 
of  soil.  In  the  top  10  mm,  the  fresh  organic  P  pool  will  be  set  to  0.0003  of 
the  initial  concentration  of  residue  on  the  soil  surface. 


orgPf^^  =  0.0003  •  rsd 


(55) 
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3  Overland  Flow  Nutrient  Processes 

Nutrient  storage  in  the  soil  occurs  primarily  in  the  uppermost  few  centi¬ 
meters  where  it  is  accessible  to  shallow-rooted  vegetation.  Surface  water 
fluxes  interact  with  this  layer  and  can  dissolve  some  of  the  soluble  nutrient 
stores,  as  well  as  cause  erosion  of  the  organic  and  adsorbed  stores.  Water 
serves  as  a  solvent  and  physical/chemical  transport  medium  for  dissolved 
and  particulate  materials  and  thus  overland  flow  can  be  an  important 
mechanism  of  nutrient  transport  from  soils  to  streams  (Sharpley  et  al. 
1992,  Goulding  et  al.  1996).  Loading  of  nutrients  from  watersheds  to  sur¬ 
face  waters  causes  changes  in  ecological  function,  and  often  has  undesir¬ 
able  environmental  and  economic  consequences.  Indeed,  non-point  runoff 
of  N  and  P  from  agricultural  landscapes  is  viewed  as  one  of  the  most 
important  factors  causing  impaired  water  quality  in  freshwater  and  estua¬ 
rine  ecosystems  (Correll  1998,  Carpenter  et  al.  1998,  Daniel  et  al.  1998, 
Downing  et  al.  1999,  Hessen  et  al.  1997,  Puckett  1995,  Sims  et  al.  1998). 
Ammonium,  nitrate  and  orthophosphate  are  the  principal  nonpoint  source 
forms  of  N  and  P.  A  major  source  for  these  elements  is  the  soil  nutrient 
pools.  Overland  flow  transport  of  these  chemicals  is  vital  for  quantifying  a 
nonpoint  source.  The  growing  concern  about  the  environmental  impact  of 
nutrients  has  enhanced  the  desire  to  predict  the  transport  and  transforma¬ 
tion  of  nitrogen  in  overland  flow  more  accurately. 

Surface  runoff  can  remove  large  quantities  of  nutrients  from  the  soil  in 
both  dissolved  and  particulate  forms.  In  the  case  of  N,  the  dissolved  inor¬ 
ganic  component  is  further  discriminated  into  NO~  and  NH'^ .  The  loss  of 

dissolved  nutrients  in  surface  runoff  is  the  result  of  rainfall  mixing  with 
the  dissolved  nutrients  in  the  upper  portion  of  the  soil.  These  dissolved 
nutrients  interact  with  surface  runoff  and  once  in  water,  they  are  trans¬ 
ported.  Suspended  nutrients,  which  are  assumed  to  be  either  organic  or 
adsorbed  inorganic  components,  attach  to  eroded  sediment  material 
derived  from  erosion  or  from  bed  erosion  in  the  stream  channel,  and  are 
transported  with  water.  The  concentration  of  suspended  nutrients  in  over¬ 
land  flow  is  largely  determined  by  soil  texture,  infiltration  rate,  rainfall 
intensity,  slope,  and  ground  cover.  However,  the  high  degree  of  temporal 
and  spatial  variability  in  the  distribution  of  surface  runoff  makes  it  diffi¬ 
cult  to  conduct  studies  on  nutrient  transport  by  this  process.  Some  of  the 
challenges  in  conducting  nutrient  mass  balance  studies  are  in  determining 
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the  unmeasured  amounts  of  nutrients  removed  from  the  water  and  stored 
in  aquatic  organisms  and  sediments,  and  nutrients  added  to  the  water 
from  resuspension,  erosion,  and  dissolution  processes. 

Nutrient  routing  in  overland  flow  is  described  by  a  two-dimensional  (2D) 
advection-dispersion  equation,  which  is  based  on  the  principle  of  conser¬ 
vation  of  mass  and  Pick’s  law.  The  forms  of  nutrient  govern  the  transport 
process  and  the  pathway  the  nutrient  is  carried  through.  The  mass  balance 
for  each  nutrient  species  in  overland  flow  in  two  directions  is  written  as: 

4  VC,  H  V  •  VC,  H  E  RjmHSj)  (56) 

12  3  4 


where 


Cj  =  concentration  of  species  j  in  overland  flow  [M/L3] 

VQ  =  5  ( •  C,  ySx  +  5  ( Qy  •  Cj  is  total  advection  transport  flux  of 


V-VQ  = 


Qx  and  Qy  = 
Dx  and  Dy= 

Rj,k  = 
Sj  = 


species  j  in  the  x-  and  y-direction  [M/L^/T] 

is  total  dispersion  flux  of  spe¬ 
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cies  j  in  the  x-  and  y-direction  [M/L^/T] 

unit  discharge  of  overland  flow  in  the  x-  and  y-direction  [L^/T] 

dispersion  coefficient  of  species  j  in  the  x-  and  y-direction 

[LVT] 

transformation  rate  of  disappearance  or  generation  of  species  j 
in  reaction  k  [M/L^/T] 

load  of  species  j  in  the  external  source  or  sink  [M/L^/T] 


In  Equation  56,  term  1  represents  advection  transport,  the  major  transport 
process  in  overland  flow.  Term  2  is  the  mass  transfer  by  diffusion  or  dis¬ 
persion.  The  third  term  is  the  sum  of  internal  mass  change  due  to  kinetic 
or  equilibrium  biological  and  chemical  transformation  processes  and 
interaction  with  the  upper  soil  layer.  All  these  processes  consider  the  effect 
on  N  and  P  cycles.  The  last  term  represents  the  external  addition  or  loss  of 
mass.  The  dissolved  substance  differs  from  the  particulate  substance  in 
overland  flow.  Within  the  model,  hydrolysis  is  defined  as  the  process  by 
which  particulate  organic  nutrients  are  converted  to  dissolved  organic 
form.  Mineralization  is  defined  as  the  process  by  which  dissolved  organic 
nutrients  are  converted  to  dissolved  inorganic  form.  Direct  mineralization 
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of  particulate  organic  nutrients  does  not  occur.  The  mass  balance  equa¬ 
tions  for  each  nutrient  variable  are  presented  below  in  detail. 

Equation  56  cannot  be  solved  alone  in  itself  but  requires  the  application  of 
a  hydrologic  model  as  an  input.  With  the  flow  field  computed  by  the 
hydrologic  model,  the  flow  results  can  be  used  to  compute  the  fate  and 
transport  of  nutrient  quantities,  i.e.,  species  concentration  and  load. 

Overland  flow  nitrogen  module 

The  dominant  N  species  in  waters  are  dissolved  inorganic  N  (DIN)  -  NH^ , 
NO^ ,  NO~ ,  dissolved  organic  N  (DON),  and  particulate  organic  N  (PON) 

and  inorganic  N  (PIN)  (Burt  et  al.  1993).  There  are  substantial  differences 
in  chemical  properties  among  the  N  species.  Total  dissolved  N  (TDN)  con¬ 
sists  of  DIN  and  DON,  and  is  readily  available  for  plant  uptake.  DIN 
mainly  comprises  NO~  +  NH'^ .  NH'^  is  the  form  of  N  taken  up  most  read¬ 
ily  by  phytoplankton  because  NO~  must  first  be  reduced  to  NH'^  before  it 

is  assimilated  into  amino  acids  in  organisms.  Organic  N  includes  all  sub¬ 
stances  in  which  N  is  bonded  to  C.  The  largest  fraction  is  made  up  of 
amino  acids  and  peptides  and  is  often  called  amino  N.  It  occurs  in  both 
particulate  and  soluble  forms.  DON  is  found  in  a  wide  range  of  complex 
chemical  forms  such  as  amino  acids,  proteins,  urea  and  humic  acids.  The 
particulate  N  consists  of  PON  and  NH'^  adsorbed  onto  mineral  particles. 

Particulate  N  can  be  found  in  suspension  or  on  the  sediment.  Some  por¬ 
tion  of  the  particulate  N  is  subject  to  rapid  mineralization,  and  is  biologi¬ 
cally  available.  Total  nitrogen  (TN)  is  a  measure  of  all  forms  of  dissolved 
and  particulate  N  present  in  waters. 

In  NSM,  N  transport  processes  are  simulated  for  dissolved  NO~ ,  dissolved 
NH'^ ,  DON  and  PON.  Figure  5  is  a  diagram  of  the  overland  flow  N  trans¬ 
port  and  transformation  module  used  in  NSM.  This  section  summarizes 
the  governing  equations  used  to  simulate  N  transport  and  transformation 
in  the  overland  flow. 
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Figure  5.  Overland  flow  nitrogen  transport  and  transformation  module. 


Particulate  organic  nitrogen 

Organic  N  may  exist  in  waters  in  particulate  form.  This  phase  includes 
small  organisms  (algae,  bacteria, ...),  both  living  and  dead,  and  fragments 
of  organisms.  As  such,  PON  particles  are  carried  by  runoff  and  transported 
through  erosion  and  deposition  processes.  The  sediment  transport  equa¬ 
tion  assumes  the  types  of  “particles”  variables  are  conservative,  which 
indicates  that  no  kinetic  functions  are  applicable.  Therefore,  mineraliza¬ 
tion,  dissolution,  or  other  transformation  processes  need  to  be  considered 
and  applied  to  PON.  Simulations  may  consider  PON  as  a  single  variable, 
or,  alternately,  represent  from  one  to  many  particle  types  or  fractions.  The 
mass  balance  equation  of  PON  transport  in  2D  overland  flow  is  written  for 
the  total  PON  concentration: 

^(h-PON  ) 

^  +  VP0/V„,  +  V .  VP0/V„,  =  Rp,,,,  - R,,,,  -  +  Spo,  (57) 


where 


PONov 

RpON,d 

RpON,e 

RpON,hyd 

SpON 


concentration  of  the  overland  flow  PON  [M/L3] 
deposition  rate  of  the  overland  flow  PON  [M/L^/T] 
erosion  rate  of  the  upper  soil  layer  organic  N  [M/L^/T] 
hydrolysis  rate  of  the  overland  flow  PON  [M/L^/T] 
external  sources  added  to  the  overland  flow  PON  [M/L^/T] 


PON  is  subsequently  entrained  into  surface  runoff  as  these  particles  are 
eroded  by  the  moving  water.  The  concentration  of  PON  available  for  trans- 
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port  is  calculated  as  the  sum  of  orgPsta,  orgPact,  and  orgPfrsh  in  the  upper 
soil  layer.  The  erosion  flux  of  PON  due  to  runoff  is  computed  by: 


R 


PON,e 


(58) 


where 

fp^  =  fraction  of  the  PON  associated  with  particle  “n” 

Vr  =  resuspension  (erosion)  velocity  [L/T] 
orgNt  =  total  concentration  of  organic  N  in  the  upper  soil  layer  [M/L3] 

When  the  sediment  transport  capacity  is  lower  than  the  sediment  load, 
sediment  deposition  occurs.  PON  particles  in  the  runoff  may  deposit  on 
the  upper  soil  layer  with  suspended  sediments.  The  deposition  flux  of  PON 
in  overland  flow  is  computed  by: 


R 


PON,d 


1 


(59) 


where 

Vse  =  effective  settling  (deposition)  velocity  [L/T] 

A  portion  of  the  PON  hydrolyzes  to  DON.  Hydrolysis  of  PON  is  modeled  as 
a  first  order  decay  process: 


^PON.hyd  ^hyd  ' 


(60) 


where 

khyd  =  hydrolysis  rate  constant  [i/T] 

Dissolved  organic  nitrogen 

While  much  research  has  studied  the  role  of  DIN  in  leachates,  the  contri¬ 
bution  of  DON  has  been  largely  ignored  or  considered  to  be  insignificant. 
However,  recent  research  demonstrated  that  DON  made  up  the  majority  of 
TDN  in  stream  exports  in  areas  with  large  anthropogenic  inputs  of  DIN. 
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The  mass  balance  equation  of  DON  transport  in  2D  overland  flow  is  writ¬ 
ten  for  the  DON  concentration: 


d{h-DON^,) 


dt 


+  VDO/V„,-V-VDO/V„, 


^PON,hyd  ^DON,f  ^DON, min  ^DON,up  ^DON 


(61) 


where 

DONov  =  concentration  of  DON  in  the  overland  flow  [M/L3] 

Rdon/  =  infiltration  rate  of  the  overland  flow  DON  [M/L^/T] 

RDON,mm  =  mineralizatiou  hydrolysis  rate  of  the  overland  flow  DON 
toN//;[M/LVT] 

Rdon,up  =  plant  uptake  rate  of  the  overland  flow  DON  [M/L^/T] 

Sdon  =  external  sources  added  to  the  overland  flow  DON  [M/L^/T] 

The  DON  infiltration  flux  can  be  computed  from  the  water  infiltration: 

RooN.f  (62) 


DON  is  mineralized  to  NH^ .  Mineralization  of  DON  is  modeled  as  a  first 
order  decay  process: 


^DON.min  ~  ^min  ' 


(63) 


where 

kmin  =  mineralization  rate  constant  [i/T] 

Ammonium  nitrogen 

The  mass  balance  equation  of  NH'^  transport  in  2D  overland  flow  is  writ¬ 
ten  for  the  dissolved  NH'^  concentration: 


dih-NH^J 


dt 


+  VNHl^-V-VNHl^ 


(64) 


=  R 


DON, min 


P  _I_P  _P  _P  _lQ 

^NHA,f  ^NH4,e  ^NHA,nit  ^NH4,up  ^NH4 


where 
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RNH4f 

RNH4,e 


RNH4,nit 

RnH4,up 

SnH4 


concentration  of  NHl  in  the  overland  flow  [M/L3] 

infiltration  rate  of  the  overland  flow  [M/L^/T] 

mass  transfer  rate  between  the  upper  soil  layer  and  surface 
runoff  [M/LVT] 

nitrification  rate  of  the  overland  flow  to  NO~  [M/L^/T] 
plant  uptake  rate  of  the  overland  flow  [M/L^/T] 
external  sources  added  to  the  overland  flow  [M/LVT] 


The  mass  transfer  flux  of  the  dissolved  between  the  soil  water  and 
the  overland  flow  can  be  expressed  as: 

^NH4,e  ~  ^ 


where 


ke  =  effective  mass  transfer  rate  constant  [L/T] 

The  dissolved  NH'^  infiltration  flux  can  be  computed  from  the  water  infil¬ 
tration: 


(66) 

Nitrification  is  represented  as  a  two-stage  process  with  NH'^  being  con¬ 
verted  to  NO^  followed  by  conversion  to  NO~ .  NO^  concentrations  are 
usually  much  less  than  NO~ ,  and  for  modeling  purposes,  NO^  is 
combined  with  NO~ .  The  nitrification  rate  is  therefore  calculated  by: 

~^nit  '^^Aov  (67) 


where 

knit  =  nitrification  rate  constant  [i/T] 
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Nitrate  nitrogen 

Nitrate  concentration  increases  with  nitrification  but  decreases  through 
denitrification.  The  mass  balance  equation  of  NO~  transport  in  2D  over¬ 
land  flow  is  written  for  the  dissolved  NO~  concentration: 


+  V/V03„,  -  V .  V/V03„, 
ot 

~  ^NH4,nit  ~  ^N03.f  ^R03,e  ~  ^N03,up  ^N03 


(68) 


where 


RN03,e 


RnosJ 

Rn03,up 

Sn03 


concentration  of  NO~  in  the  overland  flow  [M/L3] 
mass  transfer  rate  of  NO~  between  the  upper  soil  layer  and 
surface  runoff  [M/L^/T] 

infiltration  rate  of  the  overland  flow  NO~  [M/L^/T] 
plant  uptake  rate  of  the  overland  flow  NO~  [M/L^/T] 
external  sources  added  to  the  overland  flow  NO~  [M/L^/T] 


The  mass  transfer  flux  of  the  NO^  between  the  overland  flow  and  the  soil 
water  can  be  expressed  as: 

^N03,e  ~  ^  e'{^^3  ~^^3ov)  (69) 


where 


ke  =  effective  mass  transfer  rate  constant  [L/T] 

Dissolved  chemicals  in  the  runoff  will  enter  the  upper  soil  if  the  water 
transporting  those  chemicals  infiltrates.  To  account  for  this  process,  the 
NO~  infiltration  flux  can  be  computed  from  the  water  infiltration: 

RN03,=f-NO,,.  (70) 

Overland  flow  phosphorus  module 

Unlike  N,  which  is  highly  mobile,  P  solubility  is  limited  in  most  environ¬ 
ments.  P  combines  with  other  ions  to  form  a  number  of  insoluble  com- 
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pounds  that  precipitate  out  of  solution.  These  characteristics  contribute  to 
a  buildup  of  P  near  the  soil  surface  that  is  readily  available  for  transport  in 
surface  runoff.  In  addition  to  precipitating,  P  adsorbs  to  soil  particles  and 
can  be  transported  via  soil  erosion.  Sharpley  and  Syers  (1979)  observed 
that  surface  runoff  is  the  primary  mechanism  by  which  P  is  exported  from 
most  watersheds.  P  movement  in  runoff  occurs  as  particulate  P  and  dis¬ 
solved  P.  Particulate  P  is  attached  to  mineral  and  organic  sediment  as  it 
moves  with  the  runoff.  Dissolved  P  is  in  the  water  solution.  In  terms  of 
their  impact  on  eutrophication  of  water  bodies,  particulate  P  becomes  less 
available  to  algae  and  plant  uptake  than  dissolved  P  because  of  the  chemi¬ 
cal  form  it  has  with  the  mineral  (particularly  iron,  aluminum,  and  calcium) 
and  organic  compounds.  Dissolved  P  is  readily  available  for  plants,  and 
consists  of  inorganic  orthophosphate  ( H^PO~ ,  HPO^~ ,  PO^~ )  and  organic 

phosphorus-containing  compounds.  The  movement  of  dissolved  P  begins 
with  the  sorption,  dissolution,  and  extraction  of  P  from  the  soil,  plant,  and 
organic  material.  These  processes  occur  when  rain  and  runoff  water  inter¬ 
act  with  the  upper  soil  layer.  Some  water  infiltrates  into  the  soil  and  perco¬ 
lates  through  the  profile  where  desorption  of  P  will  result  in  a  low  dis¬ 
solved  concentration  in  subsurface  and  return  flow.  High  dissolved  P  con¬ 
centration  can  be  expected  in  the  water  percolating  through  organic, 
coarse-textured,  and  oxygen-depleted,  waterlogged  soils.  The  particulate  P 
consists  of  organic  P  (POP),  mineral  P,  and  phosphate  adsorbed  on  min¬ 
eral  surfaces.  Particulate  P  can  be  found  in  suspension  or  on  the  sediment. 
The  adsorption  and  desorption  of  phosphate  from  mineral  surfaces  forms 
a  buffering  mechanism  that  regulates  dissolved  phosphate  concentrations 
in  water  bodies.  Total  P  (TP)  is  a  measure  of  all  the  various  forms  of  P 
found  in  waters.  In  NSM,  P  transport  processes  are  simulated  for  DIP, 
DOP,  and  POP.  Sediment  sorbed  inorganic  P  (SIP)  will  be  added  in  future 
versions  of  NSM.  Figure  6  is  a  diagram  of  the  overland  flow  P  transport 
and  transformation  module  used  in  NSM. 

Phosphorus  is  most  commonly  assumed  to  be  transported  predominantly 
in  particulate  forms  through  soil  erosion  by  surface  runoff  (Walling  et  al. 
1997,  Bowes  et  al.  2003).  However,  where  soil  erosion  is  limited,  the 
majority  of  P  transported  by  surface  runoff  may  be  in  dissolved  forms 
(Haygarth  and  Jarvis  1997,  Heathwaite  and  Dils  2000).  Transfer  of  soil  P 
to  runoff  water  is  controlled  by  physical  and  chemical  processes  such  as 
desorption,  dissolution,  and  diffusion.  The  interaction  between  particulate 
and  dissolved  P  in  the  runoff  is  very  dynamic  and  the  mechanism  of  trans¬ 
port  is  complex.  Important  processes  related  to  P  transport  in  the  runoff 
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include:  detachment  and  deposition  of  sediment  particles;  adsorption  and 
desorption  of  soluble  P  to/from  sediment  particles  (House  et  al.  1995); 
coprecipitation  of  P  with  calcite  in  hardwaters  (House  and  Donaldson 
1986,  Jarvie  et  al.  2002);  formation  of  the  ferrous  phosphate  mineral  vivi- 
anite  in  anoxic  sediments  (Woodruff  et  al.  1999);  and  P  uptake  by  aquatic 
plants  through  either  root  or  shoot.  The  combination  of  all  of  these  proc¬ 
esses,  in  tandem  with  variations  in  overland  flow  and  other  environmental 
factors,  drives  the  following  mass  balance  equations  of  P  transport  and 
transformation  processes.  This  section  summarizes  the  governing  equa¬ 
tions  used  to  simulate  P  transport  and  transformation  in  the  overland 
flow. 


Figure  6.  Overland  flow  phosphorus  transport  and  transformation  module. 

Particulate  organic  phosphorus 

Sources  and  sinks  for  DOP  included  in  the  model  are:  hydrolysis  to  DIP, 
erosion,  deposition,  and  external  loads.  The  mass  balance  equation  of  POP 
transport  in  2D  overland  flow  is  written  for  the  POP  concentration: 

dih-POP  ) 

^ +  VPOP,„  -  V  ■  VPOP,„ 

dt  "  (71] 

_P  _p  _p  iC 

^POP,e  ^POP,d  ^POP.hyc/  ^POP 


where 

POPov  =  concentration  of  POP  in  the  overland  flow  [M/L3] 

Rpop,d  =  deposition  rate  of  the  overland  flow  POP  [M/L^/T] 

Rpop,e  =  total  erosion  rate  from  the  upper  soil  layer  organic  P  [M/L^/T] 
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Rpop,hyd  =  hydrolysis  rate  of  the  overland  flow  POP  [M/L^/T] 

Spop  =  external  sources  added  to  the  overland  flow  POP  [M/L^/T] 

Erosion  of  organic  and  adsorbed  phosphorus  occurs  in  conjunction  with 
surface  sediment  erosion  and  is  dependent  on  the  presence  of  surface  run¬ 
off.  The  concentration  of  POP  available  for  transport  is  calculated  as  the 
sum  of  orgPsta,  orgPact,  and  orgPfrsh  in  the  upper  soil.  The  erosion  flux  of 
POP  due  to  runoff  is  computed  by: 


R 


POP,e 


(72) 


where 

orgPt  =  total  concentration  of  organic  P  in  the  upper  soil  layer  [M/L3]. 
The  deposition  flux  of  POP  in  overland  flow  is  computed  by: 

N 

(73) 


A  portion  of  the  POP  hydrolyzes  to  DOP.  Hydrolysis  of  POP  is  modeled 
similarly  as: 


^POP.hyd  ^hyd  ' 


(74) 


Dissolved  organic  phosphorus 

Sources  and  sinks  for  DOP  included  in  the  model  are:  hydrolysis  from 
PON,  mineralization  to  DIP,  infiltration  into  soil,  and  external  loads.  The 
mass  balance  equation  of  DOP  transport  in  2D  overland  flow  is  written  for 
the  DOP  concentration: 


dih-DOP  ) 

— - +  VDOP„,  -  V  •  VDOP„ 

dt 

~  ^POP.hyd  ~  ^DOP.f  ~  ^DOP,min  ^DOP 


(75) 


where 

DOPov  =  concentration  of  DOP  in  the  overland  flow  [M/L3] 
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Rdop/  =  infiltration  rate  of  the  overland  flow  DOP  [M/L^/T] 

Roop.mm  =  mineralization  rate  of  the  overland  flow  DOP  [M/L^/T] 

Sdop  =  external  loadings  added  to  the  overland  flow  DOP  [M/L^/T] 

The  dissolved  organic  P  infiltration  flux  can  be  computed  from  the  water 
infiltration: 


(76) 

DOP  is  mineralized  to  phosphate.  Mineralization  of  DOP  is  modeled  simi¬ 
larly  as: 


RoOF.ivin  ~  ^min  ' 


(77) 


Dissolved  inorganic  phosphorus 

Inorganic  P  is  mainly  transported  into  the  waters  via  surface  runoff,  either 
as  particulate  P  bounded  to  sediment  particles  in  connection  with  erosion 
or  as  dissolved  P.  Once  P  is  in  solution,  the  transition  between  dissolved 
and  sorbed  forms  during  overland  flow  can  change,  mediated  by  sorption- 
desorption  properties  of  the  sediments.  A  linear  relationship  exists 
between  the  concentration  of  PIP  in  sediment  and  the  concentration  of 
DIP  in  water  at  equilibrium.  Equilibrium  partition  coefficients  are  used  to 
distribute  the  total  between  the  two  states.  The  mass  balance  equation  of 
DIP  transport  in  2D  overland  flow  is  written  for  the  DIP  concentration: 


d{h-DIPJ 

rN.  '  ’  ov  ov 

dt 

~  f^DOP.min  ~  ^DIP,f  ~  ^DIP,e  ~  ^DIP.up  ^DIP 


(78) 


where 

DIPov  =  concentration  of  DIP  in  the  overland  flow  [M/L3] 

Rdipj  =  infiltration  rate  of  the  overland  flow  DIP  [M/L^/T] 

RDip,e  =  mass  transfer  rate  between  the  upper  soil  layer  and  overland 
flow  [M/LVT] 

Rdipmp  =  plant  uptake  rate  of  the  overland  flow  DIP  [M/L3] 

Sdip  =  external  loadings  added  to  the  overland  flow  DIP  [M/L^/T] 
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The  mass  transfer  flux  of  the  dissolved  inorganic  P  between  the  soil  water 
and  the  overland  flow  can  be  expressed  as: 

(79) 

The  dissolved  inorganic  P  infiltration  flux  can  be  computed  from  the  water 
infiltration: 

(80) 
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4  In-Stream  Nutrient  Processes 

Water  quality  in  estuaries,  lakes  and  reservoirs  depends  upon  the  quantity 
and  quality  of  inflows  from  the  upstream  watershed,  which  is  usually  the 
most  significant  source  of  pollutants.  Streamflow  is  therefore  of  primary 
importance  to  the  health  of  receiving  waters  and  its  quantification  is  cru¬ 
cial  to  our  ability  to  manage  those  systems  in  an  environmentally  healthy 
manner.  In  order  to  understand  or  predict  nutrient  pollution,  the  trans¬ 
port  and  the  transformations  of  nutrients  within  the  channel  network  are 
accounted  for  in  NSM.  Parameters  that  affect  stream  water  quality  and  can 
be  considered  pollution  indicators  include  N,  P,  and  C.  The  in-stream 
nutrient  transformation  processes  formulations  that  follow  are  adapted,  in 
part,  from  the  QUAL2E  and  CE-QUAL-RIVi  (RIVi)  models.  Both 
QUAL2E  and  RIVi  are  iD  comprehensive  stream  water  quality  models.  In- 
stream  internal  sources  or  sinks  are  correlated  to  temperature,  which  gov¬ 
erns  kinetic  rates  of  biological  or  chemical  reactions  as  well  as  the  equilib¬ 
rium  speciation.  Water  quality  constituents  that  can  be  simulated  in  the 
current  NSM  version  are: 

•  Algae  as  Chlorophyll  a 

•  Organic  Nitrogen  (orgN) 

•  Ammonia  Nitrogen  ( NH'^ ) 

•  Nitrate  Nitrogen  ( NO~ ) 

•  Organic  Phosphorus  (orgP) 

•  Inorganic  Phosphorus  (irgP) 

•  Carbonaceous  Biological  Oxygen  Demand  (CBOD) 

•  Dissolved  Oxygen  (DO) 


Eigure  7  shows  the  stream  water  quality  state  variables  and  interaction. 
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Figure  7.  Schematic  of  in-stream  water  quality  state  variables  and  interaction. 


In-stream  algae 

During  the  day,  algae  increases  the  stream’s  dissolved  oxygen  concentra¬ 
tion  via  photosynthesis.  At  night,  algae  reduce  the  concentration  via  respi¬ 
ration.  As  algae  grow  and  die,  they  form  part  of  the  in-stream  nutrient 
cycle. 


Chlorophyll  a 


Chlorophyll  a  is  assumed  to  be  directly  proportional  to  the  concentration 
of  phytoplanktonic  algal  biomass. 


chla^,^a,-alg^, 


(81) 
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where 

chlach  =  chlorophyll  a  concentration  [M/V] 

Qo  =  ratio  of  chlorophyll  a  to  algal  biomass  [M/M] 
algch  =  algae  biomass  concentration  [M/V] 


Algal  growth 

Growth  and  decay  of  algae/chlorophyll  a  are  calculated  as  a  function  of  the 
growth  rate,  the  respiration  rate,  the  settling  rate,  and  the  amount  of  algae 
present  in  the  stream.  The  change  in  algal  biomass  is  given  as: 


dalggh 

dt 


Ma  '  alggh 


Pa-aiggh 


(82) 


where 

jua  =  local  specific  growth  rate  of  algae  [i/T] 

Pa  =  local  respiration  or  death  rate  of  algae  [i/T] 
cr,  =  local  settling  rate  for  algae  [L/T] 

D  =  the  depth  of  water  in  the  channel  [L] 

Local  specific  growth  rate  of  algae 

The  local  specific  growth  rate  of  algae  is  a  function  of  the  availability  of 
required  nutrients,  light,  and  temperature.  A  variety  of  mathematical 
equations  for  calculating  multiple  nutrient-light  limitations  on  algae 
growth  rate  has  been  reported.  In  QUAL2E,  the  user  has  three  options  for 
calculating  the  specific  algal  growth  rate:  1)  multiplicative;  2)  limiting 
nutrient;  and  3)  harmonic  mean. 

Multiplicative  option 

The  multiplicative  option  multiplies  the  growth  factors  for  light,  N,  and  P 
together  to  determine  their  net  effect  on  the  local  algal  growth  rate.  This 
option  has  its  biological  basis  in  the  multiplicative  effects  of  enzymatic 
processes  involved  in  photosynthesis: 


Pa, 20  Pmax 


■FL-FN-FP 


(83) 


where 

l~la,20  =  local  specific  algal  growth  rate  at  20  °C  [i/T] 
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//max  =  maximum  specific  algal  growth  rate  [i/T] 

FL  =  algal  growth  attention  factor  for  light 
FN  =  algal  growth  limitation  factor  for  N 
FP  =  algal  growth  limitation  factor  for  P 

Limiting  nutrient  option 

The  limiting  nutrient  option  calculates  the  local  algal  growth  rate  as  lim¬ 
ited  by  light  and  either  N  or  P.  The  nutrient/light  effects  are  multiplicative, 
but  the  nutrient/nutrient  effects  are  alternate.  The  algal  growth  rate  is 
controlled  by  the  nutrient  with  the  smaller  growth  limitation  factor.  This 
approach  mimics  Liebig’s  law  of  the  minimum: 

Fa.2o  mm  {FN,FP)  (84) 


Harmonic  mean  option 

The  harmonic  mean  is  mathematically  analogous  to  the  total  resistance  of 
two  resistors  in  parallel  and  can  be  considered  a  compromise  between  the 
previous  two  equations.  The  algal  growth  rate  is  controlled  by  a  multiplica¬ 
tive  relation  between  light  and  nutrients,  while  the  nutrient/nutrient  inter¬ 
actions  are  represented  by  a  harmonic  mean. 


Fa.20  /^max  ‘  ‘ 


max  * 

+ 


FN  FP 


(85) 


Aigai  growth  iimiting  factor  for  iight 

A  number  of  mathematical  relationships  between  photosynthesis  and  light 
have  been  developed.  All  relationships  show  an  increase  in  photosynthesis 
rate  with  increasing  light  intensity  up  to  a  maximum  or  saturation  value. 
The  algal  growth  limiting  factor  for  light  is  calculated  using  a  Monod  half¬ 
saturation  method.  In  this  option,  the  algal  growth  limitation  factor  for 
light  is  defined  by  a  Monod  expression: 


FL 


Z 


pht,z 


K  +/ 

^  'pht,z 


(86) 


where 

FLz  =  algal  growth  attenuation  factor  for  light  at  depth  z 
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Ipht,z  =  photosynthetically  active  light  intensity  at  a  depth  z  below  the 
water  surface  [W/L^-T],  which  is  radiation  with  a  wavelength 
between  400  and  700  nm 

Kl  =  half-saturation  coefficient  for  light  [W/L^-T],  which  is  defined 
as  the  light  intensity  at  which  the  algal  growth  rate  is 
50  percent  of  the  maximum  growth  rate 

Photosynthesis  is  assumed  to  occur  throughout  the  depth  of  the  water  col¬ 
umn.  The  variation  in  light  intensity  with  depth  is  defined  by  Beer’s  Law: 

W  =Wexp  (-/<,•  z)  (87) 


where 

Ipht,hr  =  photosynthetically  active  solar  radiation  reaching  the 

ground/water  surface  during  a  specific  hour  on  a  given  day 
[W/L2-T] 

kl  =  light  extinction  coefficient  [i/L] 
z  =  depth  from  the  water  surface  [L] 

Substituting  Equation  86  into  87  and  integrating  over  the  depth  of  flow 
gives: 


FL  = 


yk, -Dj 


In 


K  +/ 

^‘pht,hr 


exp(-k,-D) 


The  photosynthetically  active  solar  radiation  is  calculated  as: 

I  =  I  -fr 

pht,hr  hr  pht 


(88) 


(89) 


where 

Ihr  =  solar  radiation  reaching  the  ground  during  a  specific  time  on 
the  current  day  of  simulation  [W/L^-T] 
frpht  =  fraction  of  solar  radiation  that  is  photosynthetically  active 

Algal  growth  limiting  factor  for  nutrients 

The  algal  growth  limiting  factor  for  N  will  be  defined  by  a  Monod  expres¬ 
sion.  Algae  are  assumed  to  use  both  ammonia  and  nitrate  as  a  source  of 
inorganic  N. 
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FN  = 


(no,,,+nh:,,)+k. 


(90) 


where 

NO~ch  =  nitrate  concentration  in  the  stream  [M/V] 

NH'^^h  =  ammonium  concentration  in  the  stream  [M/V] 

Kn  =  Michaelis-Menton  half-saturation  constant  for  N  [M/V] 


The  algal  growth  limiting  factor  for  P  is  defined  as  a  Monod  expression. 


FP  = 


DPq, 

DPo,+K, 


(91) 


where 

DPch  =  concentration  of  dissolved  P  in  the  stream  [M/V] 

Kp  =  Michaelis-Menton  half-saturation  constant  for  P  [M/V] 

The  Michaelis-Menton  half-saturation  constant  for  N  and  P  defines  the 
concentration  of  N  and  P  at  which  algal  growth  is  limited  to  50  percent  of 
the  maximum  growth  rate.  Users  are  allowed  to  set  these  values.  Typical 
values  of  Kn  range  from  0.01  to  0.30  mg  N/L  while  K?  will  range  from 
0.001  to  0.05  mg  P/L. 

Temperature  factors  affecting  algae  growth 

Algal  growth  and  death  rates  are  temperature  dependent.  If  the  algal 
growth  rate  at  20  °C  is  calculated,  the  rate  coefficient  is  then  adjusted  for 
temperature  effects  using  a  Streeter-Phelps  type  formulation: 

(92) 


where 

Twat  =  average  water  temperature  for  the  day  [T] 

The  local  respiration  or  death  rate  of  algae  represents  the  net  effect  of 
three  processes:  1)  the  endogenous  respiration  of  algae,  2)  the  conversion 
of  algal  P  to  organic  P,  and  3)  the  conversion  of  algal  N  to  organic  N.  The 
respiration  rate  is  adjusted  to  the  local  water  temperature  using  the 
relationship: 
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A=A,2o-1.047''“-”I  (93) 

The  local  settling  rate  of  algae  represents  the  net  removal  of  algae  due  to 
settling.  The  settling  rate  of  algae  is  adjusted  to  the  local  water  tempera¬ 
ture  using  the  relationship: 

cr,  =  (7,  20  •1.024^'’™' (94) 

where 

01^20  =  local  algal  settling  rate  at  20  °C  [L/T] 

Chlorophyll  a  in  surface  runoff 

The  algal  biomass  loading  to  the  stream  can  be  estimated  as  the  chloro¬ 
phyll  a  loading  from  the  land  area.  The  chlorophyll  a  concentration  in  sur¬ 
face  runoff  is  calculated  using  a  simplified  version  of  Cluis  et  al.’s  exponen¬ 
tial  function  (1988): 

chia,,  =  0  if  (Oo,  <  10  "  m7s)  or  (tp  and  TN  <  10  ®)  (95) 

ch/a„,  =  if  (Oo,  >10  "m7s)  or  (tp  and  r/V>  10  ®)  (96) 

Qov 

=  if  (Oo,  >10  "m7s)  or  (rP<  10  ®  and  r/Vx  10  7  (97) 

Qov 

where 

chloov  =  chlorophyll  a  concentration  in  the  surface  runoff  [M/V] 

Qov  =  surface  runoff  [L] 

TN  =  total  Kjeldahl  N  load  (kmoles) 

TP  =  total  P  load  (kmoles) 

In-stream  nitrogen 

In  aerobic  water,  there  is  a  stepwise  transformation  from  organic  N  to 
ammonia,  to  nitrite,  and  finally  to  nitrate.  The  forms  of  N  simulated  by  the 
model  are  organic  N,  ammonium,  and  nitrate.  Nitrite  is  not  considered 
because  the  overall  rate  of  nitrification  is  ammonia-limited  (Parker  et  al. 
1975)?  and  stream  surveys  do  not  exhibit  significant  increases  in  nitrite  in 
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nitrifying  reaches  (Garland  1978,  Miller  and  Jennings  1979).  This  section 
summarizes  the  equations  used  to  simulate  the  N  cycle  in  the  stream. 

Organic  nitrogen 

The  amount  of  organic  N  in  the  stream  may  be  increased  by  the  conversion 
of  algal  biomass  N  to  organic  N.  Organic  N  concentration  in  the  stream 
may  be  decreased  by  the  conversion  of  organic  N  to  ammonia  or  the  set¬ 
tling  of  organic  N  with  sediment.  The  change  in  organic  N  concentration  is 
calculated  as: 


dorg/Vc,  ^  g/Vc,.o,g 
dt  dx 


=  a,-p,-a\g,. 


pN,2-orgN,,-cT^-orgN^, 


(98) 


where 

orgNch 

ai 

^N,3 

(^4 


organic  N  concentration  in  the  stream  [M/V] 
fraction  of  algal  biomass  that  is  N 
rate  constant  for  hydrolysis  of  organic  N  to  NH'^  [i/T] 
organic  N  settling  rate  [i/T] 


N^h,org  is  total  organic  N  transport  flux  in  the  channel  direction  (has 
two  components,  advective  and  dispersive): 


N 


ch,org 


=  u,-orgN^ 


-D„ 


dorgN^^ 

dx 


(99) 


Ammonium  nitrogen 

The  amount  of  ammonium  in  the  stream  may  be  increased  by  the  miner¬ 
alization  of  organic  N  and  diffusion  of  ammonium  from  the  streambed 
sediments.  The  ammonium  concentration  in  the  stream  may  be  decreased 
by  the  conversion  of  NH'^  to  NO~  or  the  uptake  of  NH'^  by  algae.  The 

change  in  ammonium  concentration  is  calculated  as: 


dN^^dr^ 

dt  dx 

=  Pn.2  •  orgN^,  -  -  Kh:  • •  a  'go. 


(100) 
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where 

(5n,i  =  nitrification  rate  constant  of  [i/T] 

0-3  =  benthos  (sediment)  source  rate  for  [M/L^-T] 

fr  =  fraction  of  algal  N  uptake  from  the  pool,  and  calculated 

as: 


f  NH^ 

'nH* 


f  NH^ 

'NHi 


■NO. 


3ch 


(101) 


where 

=  preference  factor  for  NH'^ 


Nitrate  nitrogen 

The  amount  of  nitrate  ( NO~ )  in  the  stream  may  be  increased  by  the  oxida¬ 
tion  of  NO^ .  The  conversion  of  NO^  to  NO~  occurs  more  rapidly  than 
the  conversion  of  NH'^  to  NO^ ,  so  the  amount  of  nitrite  present  in  the 

stream  is  usually  very  small  and  not  considered  in  NSM.  The  nitrate  con¬ 
centration  in  the  stream  may  be  decreased  by  the  uptake  of  NO~  by  algae. 

The  change  in  nitrate  concentration  is  calculated  as: 

I  ^^ch,N0 

dt  8x  (102) 

“  /^N,l  ■  ^^4ch  ~  /^N.2  ■  ^^3ch  ~  ~  ’  /^a  ’ 


where 

/3n,2  =  denitrification  rate  constant  [i/T] 

Temperature  factors  affecting  nitrogen  transformation 

Nitrogen  transformation  rates  are  temperature  dependent.  If  the  user 
defines  the  local  rate  constant  for  hydrolysis  of  organic  N  to  ammonia  at 
20  °C,  the  organic  N  hydrolysis  rate  will  be  adjusted  to  the  local  water  tem¬ 
perature  using  the  relationship: 


(7wa,-20) 


A/, 3  -  Pn.3, 20  •  1-047 


(103) 
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If  the  user  defines  the  rate  coefficient  for  organic  N  settling  at  20  °C,  the 
organic  N  settling  rate  will  be  adjusted  to  the  local  water  temperature 
using  the  relationship: 


0-4  =  0-4  20  •  (104) 

The  rate  constant  for  biological  oxidation  of  ammonia  N  will  vary  as  a 
function  of  in-stream  oxygen  concentration  and  temperature.  The  rate 
constant  is  calculated  as: 

A.1  =  Ai2o  ■[l-exp(-0.6-DO„)]-1.083l'--”)  (105) 


where 

Pn,i,2o  =  rate  constant  for  biological  oxidation  of  ammonia  N  at 
20  °C  (l/s) 

DOch  =  DO  concentration  in  the  stream  [M/V] 

If  the  user  defines  the  benthos  source  rate  for  ammonium  at  20  °C,  the 
benthos  source  rate  for  ammonium  N  is  adjusted  to  the  local  water  tem¬ 
perature  using  the  relationship: 

0-3  =  0-3  20  •  (106) 

The  rate  constant  for  biological  oxidation  of  nitrite  to  nitrate  varies  as  a 
function  of  in-stream  oxygen  concentration  and  temperature.  The  rate 
constant  is  calculated  as: 

=A,2,2o'[l-exp(-0.6-DO„)]-1.047"‘--“>  (107) 


where 

/3n.2,2o  =  rate  constant  for  biological  oxidation  of  nitrite  to  nitrate  at 
20  °C 

In-stream  phosphorus 

Although  P  is  not  directly  toxic  in  freshwater  systems,  it  is  often  the  limit¬ 
ing  nutrient  for  aquatic  growth  and  high  levels  can  lead  to  excessive  algal 
or  other  plant  growth.  This  can  result  in  nuisance  issues,  as  well  as  con¬ 
tribute  to  eutrophication.  Within  the  river  system,  P  undergoes  numerous 
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transformation  processes  in  the  course  of  transport.  Important  processes 
related  to  P  transformation  within  the  river  system  include,  amongst  oth¬ 
ers,  detachment  and  deposition  of  sediment  particles,  and  adsorption  and 
desorption  of  soluble  P  to/from  sediment  particles  both  in  suspension  and 
in  the  riverbed  (House  et  al.  1995). 

The  P  cycle  is  similar  to  the  N  cycle.  The  death  of  algae  transforms  algal  P 
into  organic  P.  Organic  P  is  mineralized  to  soluble  P,  which  is  available  for 
uptake  by  algae.  Organic  P  may  also  be  removed  from  the  stream  by  set¬ 
tling.  This  section  summarizes  the  equations  to  be  used  in  simulating  the  P 
cycle  in  the  stream. 

Organic  phosphorus 

The  amount  of  organic  P  in  the  stream  may  be  increased  by  the  conversion 
of  algal  biomass  to  organic  P.  Organic  P  concentration  in  the  stream  may 
be  decreased  by  the  conversion  of  organic  P  to  soluble  inorganic  P  or  the 
settling  of  organic  P  with  the  sediment.  The  change  in  organic  P  concen¬ 
tration  is  calculated  as: 

=  a,-p. .  alg„  -  .  orgP„  -  <7,  ■  orgP^,  (108) 

ot  ox 


where 

orgPch 

a2 

Ppa 


dissolved  organic  P  concentration  [M/V] 
fraction  of  algal  biomass  that  is  P  [M/M] 
rate  constant  for  mineralization  of  organic  P  [i/T] 
rate  coefficient  for  organic  P  settling  [i/T] 


Inorganic  phosphorus 

The  amount  of  soluble,  inorganic  P  in  the  stream  may  be  increased  by  the 
mineralization  of  organic  P  and  diffusion  of  inorganic  P  from  the  stream- 
bed  sediments.  The  soluble  P  concentration  in  the  stream  may  be 
decreased  by  the  uptake  of  inorganic  P  by  algae.  The  change  in  soluble  P 
concentration  is  calculated  as: 


ch,irg 


dt 


dx 


=  /3p.-orgP^,+ 


(Jn 


1000 -D 


-«2  -alg 


ch 


(109) 


where 
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irgPch  =  soluble  inorganic  P  concentration  [M/V] 

(Js  =  benthos  (sediment)  source  rate  for  soluble  P  [M/L^-T] 

Temperature  factors  affecting  phosphorus  transformation 

Phosphorus  transformation  rates  are  temperature  dependent.  If  the  user 
defines  the  local  rate  constant  for  mineralization  of  organic  P  at  20  °C,  the 
organic  P  mineralization  rate  is  adjusted  to  the  local  water  temperature 
using  the  relationship: 


44=>^p,4,20-1-047("""''°^  (110) 

If  the  user  defines  the  rate  coefficient  for  organic  P  settling  at  20  °C,  the 
organic  P  settling  rate  is  adjusted  to  the  local  water  temperature  using  the 
relationship: 


=f^5,20-l-024("“"“"°^  (111) 

If  the  user  defines  the  benthos  source  rate  for  soluble  P  at  20  °C,  the  ben¬ 
thos  source  rate  for  soluble  P  is  adjusted  to  the  local  water  temperature 
using  the  relationship: 


0-2  =  0-220  •  1.074^''™'“^°^  (112) 

In-stream  carbonaceous  biological  oxygen  demand  (CBOD) 

The  CBOD  of  the  water  is  the  amount  of  oxygen  required  to  decompose  the 
organic  material  in  the  water.  CBOD  is  added  to  the  stream  with  loadings 
from  surface  runoff  or  point  sources.  Within  the  stream,  two  processes  are 
modeled  that  impact  CBOD  levels,  both  of  which  serve  to  reduce  the  car¬ 
bonaceous  biological  oxygen  demand  as  the  water  moves  downstream.  The 
change  in  CBOD  concentration  within  the  stream  can  be  calculated  as: 

=  -r,.  CBOD,, -r,.  CBOD,,  (113) 


where 

Ki  =  CBOD  deoxygenation  rate  [i/T] 
K3  =  settling  loss  rate  of  CBOD  [i/T] 
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Temperature  factors  affecting  CBOD 

CBOD  rate  constants  are  temperature  dependent.  If  the  user  defines  the 
carbonaceous  deoxygenation  rate  at  20  °C,  the  CBOD  deoxygenation  rate 
will  be  adjusted  to  the  local  water  temperature  using  the  relationship: 

■  1.047^^"'"'“^°^  (114) 

If  the  user  defines  the  settling  loss  rate  of  CBOD  at  20  °C,  the  settling  loss 
rate  is  adjusted  to  the  local  water  temperature  using  the  relationship: 

x-3  =  x-3  20  •  1.024^''»=''^°^  (115) 


CBOD  in  surface  runoff 

The  loading  function  for  the  ultimate  CBOD  from  the  surface  runoff  is  cal¬ 
culated  based  on  a  relationship  given  by  Thomann  and  Mueller  (1987): 

CBOD,,  =0.001  (116) 

Qo,-area 


where 

CBODov  =  CBOD  concentration  in  surface  runoff  [M/V] 
orgCov  =  organic  carbon  in  surface  runoff  [M] 
area  =  area  of  a  grid  cell  [L^] 

In-stream  dissolved  oxygen 

An  adequate  dissolved  oxygen  concentration  is  a  basic  requirement  for  a 
healthy  aquatic  ecosystem.  Dissolved  oxygen  concentrations  in  streams  are 
a  function  of  atmospheric  reaeration,  photosynthesis,  plant  and  animal 
respiration,  benthic  (sediment)  demand,  biochemical  oxygen  demand, 
nitrification,  salinity,  and  temperature.  The  change  in  dissolved  oxygen 
concentration  is  calculated  as: 


tfPOc. 

dt 


=  ^2 


DO,,  )  +  {a^-g,-a^-p,)-a  Ig,,  -  ■  CBOD,, 
~  1000  D  ~  ~  ^6  ■  /^N,2  ■ 


(117) 


where 
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Kz  =  reaeration  rate  for  Fickian  diffusion  [i/T] 

«3  =  rate  of  oxygen  production  per  unit  of  algal  photosynthesis 
[M/M] 

a4  =  rate  of  oxygen  uptake  per  unit  of  algae  respired  [M/M] 

K4  =  sediment  oxygen  demand  rate  [M/L^-T] 

«5  =  rate  of  oxygen  uptake  per  unit  NHl  oxidation  [M/M] 
ae  =  rate  of  oxygen  uptake  per  unit  NO^  oxidation  [M/M] 

DOsat  =  DO  saturation  concentration  [M/V] 

An  equation  developed  by  the  American  Public  Health  Association  (APHA) 
(1985)  is  used  to  calculate  the  saturation  of  dissolved  oxygen  in 
freshwater: 


DO33,  =  exp 


-139.34410  + 


1.575701x10^  6.642308x10" 


7, 


1.243800x10 


wat,K 

10 


( '^wat.K  ) 


ij'wat.K  ) 

8.621949x10"^ 

( ~^wat,K  ) 


(118) 


where 

Tu,at,K  =  water  temperature  in  Kelvin  (273.15  +  °C) 

Reaeration 

Reaeration  occurs  by  diffusion  of  oxygen  from  the  atmosphere  into  the 
stream  and  by  the  mixing  of  water  and  air  that  occurs  during  turbulent 
flow. 

Reaeration  by  Fickian  diffusion 

Numerous  methods  have  been  developed  to  calculate  the  reaeration  rate  at 
20  °C.  A  few  of  the  methods  are  listed  below. 

Using  field  measurements,  Churchill  et  al.  (1962)  derived  the  relationship: 

x-220  =(5.03-v°-®®®-D  - - —  (119) 

V  c  j  86400 


where 


Vc  =  average  stream  velocity  [L/T] 
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O’Connor  and  Dobbins  (1958)  incorporated  stream  turbulence  character¬ 
istics  into  the  equations  they  developed.  For  streams  with  low  velocities 
and  isotropic  conditions: 
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(120) 


where 

Dm  =  molecular  diffusion  coefficient  [L^/T],  Dm  can  be  calculated: 

=177-1.037^^»“^^°^  (121) 


For  streams  with  high  velocities  and  nonisotropic  conditions: 
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where 

sip  =  slope  of  the  streambed  [L/L] 


Owens  et  al.  (1964)  developed  an  equation  to  determine  the  reaeration 
rate  for  shallow,  fast-moving  streams  where  the  stream  depth  is  0.1  to 
3.4  m  and  the  velocity  is  0.03  to  1.5  m/s. 
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Reaeration  by  turbulent  flow  over  a  dam 

Reaeration  will  occur  when  water  falls  over  a  dam,  weir,  or  other  structure 
in  the  stream.  The  amount  of  reaeration  that  occurs  is  a  function  of  the 
oxygen  deficit  above  the  structure  and  a  reaeration  coefficient: 


ADO,,  =D, 


D.=D, 


(124) 


where 

ADOch  =  change  in  dissolved  oxygen  concentration  [M/V] 
Da  =  oxygen  deficit  above  the  structure  [M/V] 
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Db  =  oxygen  deficit  below  the  structure  [M/V] 

Rea  =  reaeration  coefficient 

The  oxygen  deficit  above  the  structure  Da  is  calculated  as: 

(125) 

Butts  and  Evans  (1983)  document  the  following  relationship  that  can  be 
used  to  estimate  the  reaeration  coefficient: 

rea  =  1  +  0.38  •  coef^  •  coef,  •  •  (1  -0.11  •  )  •  (l  +  0.046  •  )  (126) 

where 

coefa  =  an  empirical  water  quality  factor 
coefb  =  an  empirical  dam  aeration  coefficient 
hfaii  =  height  through  which  the  water  falls  [L] 

Twat  =  average  water  temperature  [T] 

The  empirical  water  quality  factor  is  assigned  a  value  based  on  the  condi¬ 
tion  of  the  stream,  as  shown  in  Table  2. 


Table  2.  Empirical  water  quality  factor. 


Water  Quality  Description 

Coefa 

Clean  Water 

1.80 

Slightly  Polluted  Water 

1.60 

Moderately  Polluted  Water 

1.00 

Grossly  Polluted  Water 

0.65 

The  empirical  dam  aeration  coefficient  is  assigned  a  value  based  on  the 
type  of  structure,  as  shown  in  Table  3. 


Table  3.  Empirical  dam  aeration  coefficient. 


Structure  Type 

Coefb 

Flat  Broad  Crested  Weir 

0.70  to  0.90 

Sharp  Crested  Weir  with  Straight  Slope  Face 

1.05 

Sharp  Crested  Weir  with  Vertical  Face 

0.80 

Sluice  Gates  with  Submerged  Discharge 

0.05 
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Temperature  factors  affecting  dissolved  oxygen 

Dissolved  oxygen  rates  are  temperature  dependent.  If  the  user  defines  the 
sediment  oxygen  demand  rate  at  20  °C,  the  sediment  oxygen  demand  rate 
is  then  adjusted  to  the  local  water  temperature  using  the  relationship: 

x-4  =  x-4  20  •  1.060^''“^' (127) 

If  the  user  defines  the  reaeration  rate  at  20  °C,  the  reaeration  rate  is 
adjusted  to  the  local  water  temperature  using  the  relationship: 

x-2  =  x-220  •  1.024^''™'“^°^  (128) 


where 

x'2,20  =  reaeration  rate  at  20  °C  [i/T] 

Dissolved  oxygen  in  surface  runoff 

Rainfall  is  assumed  to  be  saturated  with  oxygen.  To  determine  the  dis¬ 
solved  oxygen  concentration  of  surface  runoff,  the  oxygen  uptake  by  the 
oxygen-demanding  substance  in  runoff  is  subtracted  from  the  saturation 
oxygen  concentration. 


DO  =D0  -CBOD  •  — 

ov  sat  1  ^  ^  ov  2^ 


(129) 


where 

DOov  =  DO  concentration  in  surface  runoff  [M/V] 

Ki  =  assumed  to  be  1.047  day  ^  for  this  formulation 

QUAL2E  and  similar  iD  models  lack  a  clear  operational  definition  of  water 
quality  parameters  within  the  model.  For  example,  it  is  known  that  many 
forms  of  organic  N  are  present  in  natural  waters.  QUAL2E  combines  them 
all  under  ‘organic  N’  and  does  not  specify  further  whether  it  is  TON, 
Kjeldahl  N,  particulate,  dissolved,  or  other.  In  contrast,  the  CE-QUAL-ICM 
(ICM)  description  contains  a  precise  specification  of  the  N,  P,  and  C  vari¬ 
ables.  ICM  is  a  finite  volume  eutrophication  model  that  can  be  used  to 
simulate  1-,  2-,  3-dimensional  water  quality  variables  including  multiple 
forms  of  algae,  carbon,  nitrogen,  phosphorus,  and  silica;  and  dissolved 
oxygen  (Cerco  and  Cole  1995).  Since  ICM  contains  more  processes,  it  also 
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contains  more  state  variables.  Future  versions  of  NSM  will  include 
additional  state  variables  in  addition  to  revised  process  descriptions  to 
better  reflect  current  understandings  of  in-stream  nutrient  processes. 
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